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ABSTRACT 


The  surf  zone  properties  and  the  on/off shore 
sediment  transport  rate  were  studied.  A  non-dimensional  surf 
zone  parameter  was  obtained  by  taking  the  ratio  of  natural 
swash  period  to  wave  period.  This  parameter*  in  addition  to 
its  wide  application  in  describing  the  breaker  types*  was 
pertinent  to  describing  the  flow  patterns  in  surf  zone  and 
defining  the  region  of  validity  of  the  similarity  solutions. 
A  similarity  model  was  applied  to  describe  the  flow  field 
of  breaking  waves  in  surf  zone.  Laboratory  results  from 
other  investigators  and  those  obtained  in  the  present  study 
were  used  to  teat  the  validity  of  the  proposed  model.  The 


results  showed  that  if  the  wave  was  of  the  spilling 
and  in  the  region  of  I,„>0.9  and  d/d^<0.6,  where  is  the 
surf  zone  parameter  defined  as  similarity 

solution  was  suitable.  A  numerical  model  was  presented  to 
predict  the  on/offshore  sediment  transport  rate.  This  model 
was  based  on  a  suspended  load  transport  equation  derived 
through  the  application  of  flow  field  and  sediment  suspen¬ 
sion  in  surf  zone.  Due  to  limited  knowledge  of  sediment 
suspension  in  surf  zone*  the  present  model  is  only  able  to 
predict  the  offshore  transport  rate.  Several  numerical 


-ii 


simulations  were  made  to  compare  with  the  laboratory  results 
of  the  author's  and  other  investigators*.  Pair  to  good 
agreement  was  observed. 
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concentration 

Fourier  constant  {'b'i/c') 

n'  c 

Fourier  constant  of  sediment  concentration 
Wave  celerity 
Drag  coefficient 

Constant  controlling  absolute  value  of  sediment 
concentration 

Wave  group  velocity 

Mean  suspended  sediment  concentration 

Maximum  suspended  sediment  concentration 

Time  function  of  suspended  sediment  concentration 

Sediment  concentration  measured  at  a  reference 
point  near  bottom 


0^(z) 


Vertical  profile  of  sediment  concentration 


xi- 


C^(z,t) 


E 


e^(t) 


g 

h 


K. 

K 


K. 


K, 


K, 


Instantaneous  sediment  concentration  in  surf  zone 

Porosity  factor 

Still  water  depth 

Median  sediment  grain  size 

Wave  energy 

Turbulence  energy  density 
Bottom  frictional  coefficient 
Gravitational  acceleration 
Local  water  depth  (d+>^  ) 

Breaking  depth 

Local  water  depth  at  i’th  grid 
Breaking  wave  height 
Incident  wave  height 
Deep  water  wave  height 


Non-dimensional  surf  zone  parameter  ( 

Proportional  constant  of  s’/ra.shing  period  (K^+Kg+K^) 
Proportional  constant  of  developing  period 
Proportional  constant  of  uprushing  period 
Proportional  constant  of  downwashing  period 


KE 


Wave  kinetic  energy 


Proportional  constant  of  kinetic  wave  energy 

Proportional  constant  of  breaking  wave  energy 

Proportional  constant  of  potential  wave  energy 

Wave  length 

Deep  water  wave  length 

Uprushing  fluid  element  mass 

Net  sediment  fliix  at  a  vertical  water  column 

Number  of  waves  in  surf  zone  (T^/T) 

Reflection  coefficient 

Wave  run-up  height 

Reynolds  number 

Slope  of  potential  energy  line 

Slope  of  total  energy  line 

Wave  period 

Swash  period 

Fluid  element  traveling  time  at  developing  stage 
Fluid  element  traveling  time  at  downwashing  stage 
Fluid  element  traveling  time  at  uprushing  stage 
Thickness  of  fluid  element 


Natural  swash  period 


-xiii 
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V  Swash  velocity  of  fluid  element 

V’  Non-dimensional  velocity  (V/V* ) 

Vo  Initial  uprush  velocity  at  SWL 

Vq  Fluid  element  velocity  at  developing  stage 

V^  Fluid  element  velocity  at  downwash  stage 

Fluid  element  velocity  at  uprush  stage 

^d(max)  do^vnwash  velocity  at  SWL 

X  Fluid  element  traveling  distance 

X'  Non-dimensional  distance  (X/Xo) 

X^  Non-frictional  traveling  distance  (V  /ZgsinoO 

X  „  Maximixm  uprush  distance 

max  ^ 

z  Reference  coordinate  of  sediment  concentration 

®  measurement 

OL  Beach  slope  angle 


“Xiv 


r 

r< 

r 

fs 

CO 

K 

n 

n 

O' 

6 

T 

m 


Mean  of  wave  height  to  water  depth  ratio 

Battjes'  surf  zone  parameter 

Incipient  breaking  criteria 

Galvin's  inshore  parameter 

Fluid  density 

Sediment  density 

Particle  falling  velocity 

Breaking  wave  height  to  water  depth  ratio 

Water  elevation 

Wave  set-up 

Fundamental  wave  frequency  (2  /T) 

Kinematic  eddy  viscosity 

Bottom  shear  stress 

Shear  velocity 

Void  ratio 

Beach  slope 

n^^  wave  phase  angle 

n^^  velocity  phase  angle 

n^^  sediment  concentration  phase  angle 


CHAPTER  1 
INTRODUCTION 


The  surf  zone  properties  and  the  on/offshore  sedi¬ 
ment  transport  process  were  closely  examined  in  this  study 
to  provide  both  theoretical  and  numerical  predictions  of 
breaking  wave  properties  and  sediment  transport  rate.  The 
surf  zone  properties  studied  included  the  surf  zone  para¬ 
meter  and  the  breaking  wave  similarity.  These  properties 
were  developed  theoretically  and  calibr*ated  with  laboratory 
data.  The  laboratory  experiments  conducted  included  water 
elevation  and  mean  horizontal  velocity  measurements  of 
breaking  waves  and  the  instantemeous  suspended  sediment 
concentration  measurements  in  surf  zone.  A  numerical  model 
was  then  developed  based  on  the  theoretical  background  to 
predict  the  on/offshore  sediment  transport  rate. 

In  Chapter  2  a  comprehensive  review  of  all  related 
previous  works  were  presented?  these  included  surf  zone 
properties!  on/offshore  sediment  trsmsport  process  and  a 
review  of  laboratory  measuring  techniques. 
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The  theoretical  background  was  introduced  in  Chapter 
3.  A  non-dimensional  surf  zone  parameter  was  developed  first 
based  on  the  swash  mechanics  of  a  single  breaking  wave.  The 
physical  implications  of  this  parameter  and  its  role  in 
characterizing  the  breaking  wave  properties  were  discussed 
in  detail.  The  breaking  wave  similarity,  including  the 
similar  properties  of  water  elevation  and  horizontal  velo¬ 
city,  were  developed  based  on  the  concept  of  harmonic 
transfer  in  shoaling  waves.  A  sediment  transport  equation 
was  developed  based  on  the  derived  velocity  field  and  empi¬ 
rical  suspended  sediment  concentration.  This  transport 
equation  was  then  applied  to  a  numerical  model  for  on/off¬ 
shore  sediment  transport  rate  prediction. 

The  procedures  and  arrangements  of  laboratory 
experiments  are  discussed  in  Chapter  4.  All  tests  were 
conducted  in  the  wave  tank  located  in  the  Civil  Engineering 
Laboratory.  Both  pleuie  and  sand  beaches  were  applied.  There 
were  2  sets  of  tests  conducted,  including  10  runs  of  break¬ 
ing  wave  similarity  measurement  and  4  runs  of  suspended 
sediment  measurement.  The  air  bubble  releasing  technique  and 
photographic  method  were  applied  for  velocity  measurement 
and  the  Iowa  Sediment  Concentration  Measuring  System  was 
used  for  suspended  sediment  measurement. 


-3- 


In  Chapter  5  the  results  of  analyzed  laboratory  data 
and  the  res\ilts  of  numerical  prediction  were  presented.  The 
technique  of  Fourier  Analysis  was  applied  to  analyze  the 
harmonic  components  of  vrater  elevation  and  horizontal 
velocity.  The  author's  resvilts  were  compared  with  the  data 
from  previous  investigators  and  the  region  of  applicability 
of  similarity  property  were  discussed.  Spectral  analysis  and 
data  averaging  techniques  were  employed  for  suspended  sedi¬ 
ment  study.  The  correlation  of  sediment  concentration  and 
breaking  wave  was  studied  based  on  the  analyzed  results  of 
suspended  sediment  concentration  and  corresponding  breaking 
wave  property.  Finally  the  simulated  numerical  results  of 
the  on/offshore  sediment  transport  model  were  compared  with 
the  author's  and  the  previous  researchers'  laboratory  data. 
Sensitivity  analysis  and  prototype  simulation  of  the 
numerical  model  were  also  conducted. 

Conclusions  were  presented  in  Chapter  6.  It 
concluded  the  findings  and  the  physical  importance  of 
theoretical,  laboratory  and  numerical  res\ilts  of  the  surf 
zone  properties  and  the  on/offshore  sediment  transport 
process. 


The  measured  water  elevation  and  horizontal  velocity 
data  were  presented  in  Appendix  A  in  non-dimensional  form. 

A  theoretical  analysis  of  the  dynamics  of  air  bubble  in 
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unsteady  wave  field  was  presented  in  Appendix  B  to  demon¬ 
strate  the  percentage  error  of  the  present  velocity  measur¬ 
ing  technique.  Appendix  C  summarized  the  field  and 
laboratory  results  of  suspended  sediment  concentration  in 
surf  zone.  Both  theoretical  prediction  and  laboratory  wave 
set-ups  were  shown  in  Appendix  D.  A  list  of  computer  program 
was  presented  in  Appendix  E. 


CHAPTER  2 
LITERATURE  REVIEW 


The  previous  research  works  in  the  fields  of  surf 
zone  properties)  on/offshore  sediment  transport,  and 
laboratory  measuring  techniques  were  reviewed.  The  resear¬ 
chers  who  made  important  contributions  were  outlined  and 
their  findings  were  summarized.  Among  the  research  results, 
many  empirical  and  experimental  conclusions  were  found  to  be 
valuable  towards  the  theoretical  development  of  the  present 
study.  The  following  sections  discuss  the  literature  review 
in  detail. 

2.1  Surf  Zone  Properties 

Due  to  the  complicated  natxire  of  breaking  wave 
induced  turbulence  in  surf  zone,  few  Euialytical  solutions 
can  be  applied  successfully  in  this  region.  However,  many 
previous  empirical  formula  and  laboratory  results  indicated 
that  there  exist  some  important  surf  zone  properties  which 
woxild  lead  to  a  better  understanding  of  the  flow  field  in 
surf  zone.  The  major  surf  zone  properties  include  the  flow 
pattern,  the  surf  zone  parameter,  and  the  breaking  wave 
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slmilarity.  These  properties  are  discussed  in  detail  in 
the  following  sub-sections. 

2.1.1  Flow  Pattern  in  Surf  Zone 

The  fluid  flow  in  surf  zone  can  be  categorized  into 
three  conditions  based  on  the  degree  of  interference  among 
breaking  waves.  Following  Kemp  (1968),  the  wave-beach 
relationship  can  be  classified  into  three  categories! 

(1)  Swash  Conditioni  The  broken  wave  is  able  to 
surge  up  to  run  up  limit  on  beach  and  return  as 
backwash  to  the  breaking  point  before  the  succ¬ 
eeding  wave  had  broken.  The  flow  shoreward  the 
breaker  is  distinctly  oscillatory. 

(2)  Transition  Condition!  The  backwash  wave  can  not 
be  completed  before  the  next  wave  breaks.  The 
condition  involved  interaction  between  the  back¬ 
wash  wave  and  the  next  uprush  wave,  resulting  in 
the  interchange  of  water  between  the  zones  sea¬ 
ward  and  shoreward  of  breakers. 

(3)  Surf  Conditions  The  transition  flow  give  way  to 
flow  conditions  in  which  successive  lines  of 
breakers  continuously  spills  water  into  the 
surf  zone,  producing  a  seaward  return  flow. 


The  sequence  of  flow  conditions  is  illustrated  in 
Figure  2-1. 
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Figure  2>1  Breaker  Classification  Based  on  Surf  Zone 
Parameter  {  ) 
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The  degree  of  interference  among  breaking  waves  is 
also  closely  related  to  the  flow  pattern.  Consider  waves 
breaking  on  a  plane  beach  with  constant  slope.  During  swash 
condition,  the  broken  wave  would  complete  its  journey  up  to 
the  run  up  limit  and  return  to  breaking  point  before  the 
next  wave  breaks.  There  is  no  interference  among  breaking 
waves.  The  mean  seaward  return  flow  is  zero  in  this  con¬ 
dition. 

As  the  period  of  incident  wave  gets  shorter,  the 
number  of  breaking  waves  passing  through  a  control  section 
would  be  more  than  unity.  The  interference  between  con¬ 
secutive  breaking  vraves  results  in  an  interchange  of  fluid 
mass  seaward  and  shoreward  of  the  control  sections  in  surf 
zone. 

2.1.2  Surf  Zone  Parameters 

In  the  surf  zone,  a  non-dimensional  parameter, 
a  combination  of  beach  slope  and  incoming  wave  steepness 
has  been  shown  to  be  a  governing  factor  of  characterizing 
breaking  wave  properties.  A  review  of  the  previous  related 
works  was  made.  The  investigators  who  developed  and  applied 
the  surf  zone  parameter  included  Iribarren  and  Nogales, 

Hunt,  Bowen  et.  al. ,  Galvin,  and  Battjes,  The  various  forms 
of  the  surf  zone  parameter  proposed,  how  they  were  arrived 
at,  and  the  breaking  wave  properties  characterized  by  this 
parameter  are  summarized. 
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Iribarren  and  Nogales  (19^9)  first  used  this 
parameter  to  determine  whether  wave  breaking  would  occur. 

By  applying  the  shallow  water  trochoidal  theory  for  uniform, 
progressive  waves,  the  condition  of  incipient  breaking  can 


be  expressed  ast 


(i-i) 


where  is  the  incident  wave  height,  d^  is  the  still  water 
depth,  and  "c"  denotes  the  incipient  breaking.  The  depth  at 


breaking  was  expressed  as: 

-  ‘a"  ^  OCc 


(2-2) 


where  is  the  breaking  wave  length  expressed  as  1^,=  T^^gd^' 
and  tan^  is  the  beach  slope,  is  the  wave  period. 


Therefore, 


C2-3) 


Elimination  of  d  between  (2-1)  and  (2-3)  gives 

( “7" 

Substituting  in  (2-b)  where  is  the  deep  water 

wave  length  and  rearranging  gives 
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The  derivation  given  by  Iribarren  and  Nogales  sug¬ 
gested  that  the  incipient  breaking  of  periodic  waves 
correspond  to  a  critical  value  of  the  proposed  surf  zone 
parameter  defined  in  Equation  (2-5). 

Hmt  (1959)  suggested  that  the  wave  run-up  height 
is  proportional  to  the  surf  zone  parameter.  By  defining  the 
rvin-up  height  "R"  as  the  maximum  elevation  of  the  waterline 
above  the  undisturbed  water  level,  he  proposed  the  empirical 
theory  for  the  run-up  height  of  waves  breaking  on  a  smooth 
slope  asj 

where  is  a  porosity  factor,  approximately  1  for  solid 
beaches,  Ho  is  the  deep  water  wave  height,  and  L#  is  the 
deep  water  wave  length.  This  empirical  formula  was  proved 
to  be  in  good  agreement  with  the  observed  laboratory  results 
by  Bowen  et.  al.  (1968). 

In  addition  to  wave  run-up,  Bowen  et.  al.  (1968) 
further  pointed  out  that  the  maximum  set-up  is  also 
governed  by  the  surf  zone  parameter.  The  experimental 
results  summarized  by  them  indicated  that  the  mean  of  the 
ratio  of  wave  height  to  water  depth  in  the  surf  zone  (  /*  ) 
is  a  function  of  the  surf  zone  parameter.  Since  the  coeffi¬ 
cient  of  proportionality  of  wave  set-up  is  a  function  of  |r  , 
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the  maximiim  height  of  wave  set-up  is  governed  by  the  surf 
zone  parameter. 


Galvin  (1968,  1972)  proposed  that  the  transition 
ranges  of  the  main  breaker  types  are  also  characterized  by 
the  surf  zone  parameter.  By  introducing  the  "Offshore 
Parameter" ,  ,  the  transition  values  separating 

the  main  breaker  types  were  givens 


Surging  or  Collapsing  if 

Plunging  if  <4iS  (2-j) 

Spilling  if 

The  "Inshore  Parameter" ,  ,  presented 

by  Galvin  is  not  equivalent  to  the  surf  zone  parameter. 
However,  a  re-analysis  of  Galvin’s  data  in  terms  of 

by  Battjes  (197^)  showed  that  the  classi¬ 
fication  of  breaker  types  could  be  performed  equally  well 
with  as  with  Galvin's  inshore  parameter.  The  following 

approximate  transition  values  were  founds 


Surging  or  Collapsing 

if 

Plunging 

if 

(2-8) 

Spilling 

if 

Battjes  (197^)  re-derived  this  surf  zone  parameter 


through  dimensional  analysis  and  summarized  all  breaking 
wave  properties  characterized  by  this  parameter.  The 
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derivation  was  based  on  the  flow  parameters  of  a  periodic, 
long-crested  wave  approaching  a  rigid,  plane,  impermeable 
beach.  The  motion  of  this  periodic  wave  on  beach  was  deter¬ 
mined  by  a  dimensionless  dependent  variable  as: 

%  =  «-w 

where  is  the  beach  slope  angle,  H  is  the  incident  wave 
height  at  the  toe  of  the  beach  slope,  d  is  the  still  water 
depth,  is  the  deep  water  wave  length  and  Re  is  a  typical 
Reynolds  Number. 

For  breaking  waves,  the  values  of  Reynolds  Number 
and  relative  depth  would  not  significantly  affect  the 
resultant  motion.  Therefore,  he  concluded  that  for  many 
overall  properties  of  breaking  waves  the  dependent  variables 
were  reduced  to: 

?c~  oe 


The  characteristic  values  of  the  surf  zone  parameter  corres¬ 
ponding  to  the  overall  properties  of  breaking  waves  in  surf 
zone  was  summarized  by  him  as  shown  in  Table  2-1.  Also, 
Battjes  suggested  that  5"**  approximately  proportional 
to  the  number  of  wave  lengths  within  surf  zone. 

In  svimmary,  the  previously  developed  surf  zone 
parameters  were  arrived  empirically  or  through  physical 


Table  2-1  Summary  of  Breaking  Wave  Properties  Based  on  Surf  Zone  Parameter 
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modeling.  Although  presented  in  different  forms i  this 
parameter  can  he  generalized  as  a  ratio  of  the  beach  slope 
to  the  square  root  of  the  incident  wave  steepness.  The 
breaking  wave  properties  proved  to  be  characterized  by  the 
surf  zone  parameter  included  the  breaking  wave  criterion, 
the  breaker  type,  the  relative  importance  of  wave  run-up 
and  set-up,  the  breaker  height- to-depth  ratio,  the  number 
of  waves  in  surf  zone,  and  the  reflection  coefficient  on 
beach. 

2.1.3  Breaking  Wave  Similarity 

When  a  wave  moves  into  shoaling  water,  it  will  grow 
in  wave  height  imtil  the  water  depth  is  shallow  enough  to 
initiation  of  breaking.  The  assumption  made  for  a  shoaling 
wave  is  that  there  is  no  loss  of  energy  and  no  reflection. 

It  is  assumed  that  the  wave  period  remains  constant  in  water 
of  any  depth,  whereas  the  wave  length,  velocity  and  height 
vary,  A  realistic  wave  form  is  better  described  by  non¬ 
linear  wave  theories.  For  a  single  frequency  laboratory  wave 
shoaling  on  beach,  its  energy  would  be  transformed  from 
primary  to  higher  harmonics. 

It  was  shown  by  Flick  (1978)  that  the  transfer  is 
well  described  by  third  order  Stokes  theory  that  conserves 
energy  flux  as  long  as  the  local  Ursell  Number  remains  less 
than  1.  Depending  upon  the  bi*eaker  type,  more  or  less  energy 
relative  to  the  primary  is  transferred  to  the  higher 
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harmonics  in  the  shoaling  process  before  breaking.  The 
amplitude  dispersion  among  harmonics  together  with  phase 
lag  due  to  differential  shoaling  of  each  harmonic  component 
result  in  a  highly  asymmetric  wave  form  upon  wave  breaking. 

After  wave  breaking,  the  broken  wave  at  seaward 
limit  of  surf  zone  remain  asymmetric.  This  asymmetric 
property  is  a  result  of  relative  importance  of  higher  har¬ 
monic  amplitude  and  the  introduction  of  harmonic  phase  due 
to  wave  breaking.  The  shape  of  a  breaking  wave  is  determined 
by  the  relative  hazmonic  amplitudes  and  harmonic  phases. 

Flick  (1978)  measured  the  breaking  wave  height  and 
made  the  harmonic  analysis  as  well.  It  was  shown  that  the 
harmonic  amplitudes  as  well  as  the  spilling  wave  heights 
decay  at  a  constant  rate  with  local  water  depth  as  shown  in 
Figure  2-2.  The  result  of  harmonic  analysis  by  Flick  implies 
that  the  non-dimensional  ratio  of  breaking  wave  height  to 
local  water  depth  is  bearing  a  similar  shape. 

During  shoaling  process,  the  transfer  of  wave  energy 
to  higher  harmonic  components  also  introduces  an  asymmetric 
velocity  profile.  The  corresponding  asymmetric  velocity 
profile  was  evidenced  by  Iwagaki  (1970)  and  Flick  (1978)  as 
shown  in  Figures  2-3  and  2-4. 

The  bottom  friction  during  breaking  process  plays 
an  important  role  in  energy  dissipation.  For  higher  harmonic 


Note;  non-dimensional  horizontal  velocity 

1/H  !  non-dimensional  water  elevation 


Figure  2-3  Asymmetric  Water  Elevation  and  Horizontal 
Velocity  Profiles  by  Iwagaki  (1970) 
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Figure  2-4  Asymmetric  Water  Elevation  and  Horizontal 
Velocity  Profiles  ty  Flick  (1978) 


♦t 


-19- 


velocity  components  with  high  oscillating  frequency,  the 
frictional  dissipation  is  greater,  resulting  in  less  energy 
content  in  the  higher  harmonics.  The  smaller  contribution 
of  higher  hannonic  components  results  in  a  relatively  more 
symmetric  velocity  profile  than  the  corresponding  wave 
profile , 

Laboratory  and  field  observations  by  Flick  (1978) 
and  Thornton  ( 1976)  both  suggested  that  the  shape  of  the 
horizontal  velocity  profile  is  more  symmetric  than  the  shape 
of  the  corresponding  wave  profile.  Also,  there  exists  a  time 
lag  between  the  velocity  and  wave  profile. 

2.2  On/offshore  Sediment  Transport 

Sediment  transport  is  the  movement  of  sediment  due 
to  the  interactions  of  sediment  with  wave,  current,  tide, 
wind  and  other  causative  forces  in  the  surf  zone.  The 
direction  of  the  sediment  movement  can  be  alongshore  or 
on/offshore,  which  are  known  as  longshore  transport  auid 
on/off shore  sediment  transport,  respectively. 

The  knowledge  of  on/offshore  sediment  transport 
criteria,  transport  mechanisms,  and  modes,  and  the  response 
of  various  beach  materials  is  essential  for  the  under¬ 
standing  of  transport  process.  Previous  laboratory  experi¬ 
ments,  field  studies,  and  numerical  models  have  provided 
some  qualitative  as  well  as  qiiantitative  results  and 
empirical  formula  in  describing  the  on/offshore  transport 
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process.  These  works  are  discussed  in  the  following 
sub-sections. 

2.2.1  On/offshore  Transport  Criteria 

It  is  well  known  that  the  seasonal  variations  cf 
beach  profiles  are  associated  with  seasonal  variations  in 
wave  climate.  Specifically,  the  idealized  winter  or  sto'r 
profile  is  characterized  by  a  mild  foreshore  slope  and  the 
presence  of  an  offshore  sand  bar.  The  idealized  summer  or 
normal  profile  is  of  steeper  foreshore  slope  and  no  offshore 
bar  is  present.  These  typical  summer  and  winter  profiles 
are  illustrated  in  Figiire  2-5.  The  transition  of  summer  and 
winter  profiles  indicates  the  transition  of  on/offshore 
transport  direction. 

Laboratory  measurements  indicated  that  the  wave 
steepness  transition  between  summer  and  winter  profiles  is 
approximately  0.025.  Field  observations  also  suggested  that 
the  variation  of  profiles  are  due  to  the  seasonal  changes  in 
deep  water  wave  steepness  defined  as  H,/L,  .  Saville  (1957) 
noted  that  for  large  scale  tests,  steepness  smaller  than 
0.025  may  result  in  a  storm  profile  too.  Iwagaki  and  Noda 
(1962)  found  that  the  ratio  of  wave  height  to  sand  grain 
size  was  also  an  important  parameter.  Dean  (1973)  related 
the  on/offshore  transport  criteria  to  both  the  deep  water 
wave  steepness  and  the  sediment  falling  velocity. 


Figure  2-5  Typical  Summer  and  Winter  Profiles  and  the  Definition 
of  Energy  Regions  on  Beach 
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In  his  transport  model,  Dean  (1973)  assumed  that 
the  action  of  breaking  waves  is  sufficient  to  place  sand 
into  suspension  over  at  least  a  portion  of  the  water  column. 
After  the  suspension  phase,  and  during  the  return  of  the 
sediment  particle  to  the  bottom,  the  direction  of  net 
horizontal  displacement  depends  on  whether  the  sand  particle 
is  acted  upon  predominantly  by  onshore  or  offshore  fluid 
particle  velocities.  If,  on  the  other  hand,  the  fall  time 
is  long  compared  to  the  wave  period,  the  predominant  sedi¬ 
ment  particle  is  directed  offshore.  The  classification  of 
normal  and  storm  profiles  proposed  by  Dean  is  shown  in 
Figure  2-6.  In  an  attempt  to  relate  the  on/offshore  trans¬ 
port  criteria  to  the  surf  zone  parameter  the  author  also 
compiled  a  similar  normal  and  storm  profiles  classification 
curve  as  shown  in  Figure  2-7.  It  can  be  seen  that  the  surf 
zone  parameter  also  plays  an  important  role  in  character¬ 
izing  the  on/offshore  transport  criteria. 

2.2.2  Modes  and  Mechanism  of  Transport 

Sediment  on  beach  can  be  moved  on/off shore  by  bottom 
friction  due  to  wave  induced  bottom  fluid  particle  velocity 
and/or  on/offshore  currents.  Once  the  sediment  is  loosened 
by  fluid  particle,  it  can  be  transported  on  or  offshore  in 
the  dominant  force  direction.  This  mode  of  sediment  movement 
is  usually  called  bed  load  transport. 

Eagleson  et.  al.  (1963)  proposed  a  "null  point" 
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Figure  2-7  Classification  of  Normal  and  Storm  beach  profiles 
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analysis  which  involved  the  estimation  of  lift,  drag  and 
inertial  forces  on  uniform  sand  size  of  spherical  shape. 

The  null  point  represents  a  balance  point  between  wave 
forces  and  gravity.  Shorewards  of  this  point  the  sediments 
move  onshore  whereas  seaward  of  it  sediment  moves  offshore. 
This  null  point  analysis  renders  no  practical  value,  since 
the  particles  are  actually  irregular  in  shape. 

Another  mechanism  of  wave  induced  sediment  motion 
can  be  explained  by  the  asymmetric  bottom  shear  stress  due 
to  the  asymmetric  fluid  velocity  field.  The  net  on/offshore 
thrust  would  result  in  a  net  shoreward/seaward  movement  of 
sediments. 

Sediments  can  also  be  carried  on/offshore  by 
currents  as  suspended  load.  This  mode  of  transport  includes 
the  dynamic  processes  of  sediment  entrainment  and  trans¬ 
portation.  The  entrainment  process  requires  force-producing 
near  bottom  fluid  velocities  and  accelerations  in  excess  of 
some  threshold  value,  while  transportation  processes  require 
the  sediment  particle  moves  with  the  fluid.  This  process  can 
be  approximated  where  the  sediments  are  small  enough  to  be 
in  the  Stokes  range  of  settling  with  settling  velocity 
Reynolds  Number  less  than  1. 

The  major  sediment  carrying  mechanism  of  suspended 
load  transport  is  nearshore  current.  It  was  pointed  out  by 
Longuet -» Higgins  (1953)  Ippen  (1966)  that  when  waves 


are  approaching  a  straight  and  parallel  shoreline,  the 
continuity  of  mass  necessitates  a  vertical  distribution  of 
mass  transport  foreward  at  the  surface  and  near  bottom, 
return  flow  near  mid-depth  before  wave  breaking.  There  were 
considerable  laboratory  evidences  supporting  the  two-dimen¬ 
sional  nearshore  circulation  model  as  illustrated  in 
Figure  2-8(a).  In  breaker  zone,  the  two-dimensional  flow 
field  can  be  described  as  oscillatory,  transition  and 
circulatory  conditions.  During  circulatory  flow  condition 
the  growing  strength  of  near  bottom  return  flow  combined 
with  sediment  suspension  due  to  breaking  turbulence  results 
in  suspended  load  in  offshore  direction. 

In  the  field,  the  near-shore  current  field  is  never 
restricted  in  two-dimensional.  A  longshore  variation  of 
energy  density  and/or  an  irregular  nearshore  bottom  topo¬ 
graphy  would  induce  high-energy  concentrated  rip  currents. 

Bowen  (1969)  developed  a  theoretical  model  using  the 
concept  of  radiation  stress.  A  three-dimensional  nearshore 
circulation  model  was  proposed  by  Shepard  and  Inman  (1951) 
as  shown  in  Figure  2-8(b).  It  was  emphasized  that  the  rips 
play  a  major  role  in  carrying  the  suspended  sediment  further 
out  of  the  offshore  submarine  bar. 

2.2.3  Beach  Materials 

Sediment  transport  is  the  combined  result  of  flow 
field  in  the  nearshore  region  and  the  dynamic  response  of 
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Figure  ZS  Two-  and  Three-dimensional  Nearshore 
Circulation  Model  After  Ippen(1966) 


-28- 


beach  materials,  therefore,  the  physical  characteristics  of 
beach  materials  also  play  an  important  role  in  the  sediment 
transport  process. 

The  properties  of  beach  materials  include  their 
size,  density,  shape,  and  permeability.  Large  materials  such 
as  shingle  and  gravel  is  gradually  rolled  along  the  bottom. 
Finer  material  such  as  sand  can  be  moved  by  suspension.  The 
physical  characteristics  of  sand  important  to  its  transport 
process  are  the  grain  size  distribution,  the  submerged 
density  in  sea  water,  and  the  sphericity  and  roundness  of 
grains.  These  three  properties  can  be  combined  into  one 
important  physical  parameter  known  as  the  particle  falling 
velocity  . 

Dean  (1973)  suggested  that  this  parameter  is  vital 
for  on/offshore  transport  criteria.  Sunamura  and  Horikawa 
(197^)  indicated  that  the  ratio  of  sand  size  to  wave  length 
is  an  important  physical  parameter  controlling  sand  beach 
transportation.  Under  the  same  deep  water  wave  condition 
and  initial  beach  slope,  their  laboratory  results  showed 
erosion  for  beach  with  fine  grain  size  and  accretion  for 
beach  with  coarse  grain  size. 

2.2.4  Numerical  Modeling 

The  knowledge  of  on/offshore  transport  rate  as  well 
as  the  beach  profile  evolution  is  essential  for  predicting 
the  short  terra  beach  changes  such  as  daily  variation  of 
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beach  profile,  the  response  under  extreme  wave  condition, 
and  the  time  needed  to  recover  the  eroded  foreshore  berm 
crest.  It  is  also  useful  to  other  coastal  engineering 
problems  including  beach  nourishment,  seawall  design,  etc. 

In  an  attempt  to  model  the  on/offshore  sediment 
transport  process,  the  model  developed  should  have  the 
following  capabilities  to  be  qualified  as  "good"  models 
(Dally  and  Dean,  1980)  s 

(1)  Generate  profiles  of  both  the  normal  and  storm 
types  depending  on  the  wave  conditions  and 
sediment  characteristics. 

(2)  Predict  the  proper  shape  of  these  profiles* 
i.e.,  (a)  the  normal  profile  should  be  monotonic 
and  concave  upwards  and  (b)  the  bar  of  the 
storm  profile  should  have  the  proper  spacing 
and  shape. 

(3)  Correctly  predict  the  rate  of  profile  evolution. 

(4)  Respond  to  changes  in  water  level  due  to  tides, 
storm  surges,  or  long-term  fluctuations. 

(5)  Approach  an  equilibrium  if  all  the  relevant 
parameters  are  held  constant. 

There  has  been  plenty  of  information  published  on 
beach  profiles,  and  a  great  number  of  studies  undertaken  to 
predict  the  condition  for  bar  formation  and  profile  erosion 
or  accretion.  However,  there  is  still  a  lack  of  a  good  model 
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for  prediction  of  profile  evolution  and  rates  of  shoreline 
erosion,  although  there  have  been  a  few  good  models  that 
were  developed  for  limited  use  which  require  both  cali¬ 
bration  for  each  use  and  large  data  base  for  statistical 
analysis. 

Fox  and  Davis  (1973)  assumed  an  arbitrary  form  for 
the  beach  profile  and  manipulated  its  coefficients  empi¬ 
rically  to  change  the  bottom  with  time.  The  model  requires 
substantial  atmospheric  data  to  calibrate  and  is  extremely 
site  specific,  although  they  achieved  good  results.  Swart 
( 197^)  developed  a  transport  model  which  represents  the 
sediment  movement  from  the  inshore  region  to  the  offshore 
region  in  a  gross  sense.  The  beach  is  divided  into  two 
horizontal  layers  with  continuity  being  satisfied  between 
the  two.  Suspended  sediment  transport  is  the  dominant 
mechanism  and  the  model  can  be  adapted  to  various  wave  and 
sediment  characteristics.  Wang,  Dalrymple  and  Shiau  (1975) 
modeled  sediment  transport  by  dividing  the  nearshore  region 
at  the  breaker  line.  Inside  the  surf  zone,  a  gross  longshore 
transport  is  proportioned  across  the  region  with  the  shape 
of  the  longshore  current  determined  by  Longuet-Higgins. 
Outside  the  surf  zone,  transport  is  determined  by  super¬ 
imposing  an  exponential  sediment  concentration  profile  onto 
a  mean  flow  profile,  which  gives  good  results  for  regions 
far  from  the  breaker  line. 
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Fleming  and  Hiint  (1976),  using  both  bed  load  and 
suspended  load  mechanics,  formulated  a  three-dimensional 
transport  model  for  a  specific  beach  segment.  Profile 
changes  were  caused  by  deviations  in  the  longshore  current 
much  like  the  model  by  Wang  et.  al. ,  (1975)*  Kana(1979) 
measured  sediment  concentration  under  fully  plunging 
breakers  on  a  stable  monotonic  beach.  In  the  model,  the 
nearshore  region  is  artificially  divided  into  sections  where 
the  transport  is  governed  by  different  equations.  The  trans¬ 
port  inside  surf  zone  is  coupled  to  the  transport  calculated 
just  outside  the  breaker  line  using  Wang’s  model  which  is 
valid  only  in  regions  far  from  the  surf  zone.  The  model 
does  approach  a  state  of  equilibrium  amd  is  applicable  to 
various  sediments  and  wave  characteristics  which  makes  it 
attractive. 

Winant,  et.  al,  (1975)  described  the  variability  in 
a  beach  profile  by  means  of  empirical  eigenfunctions  with 
good  results.  However,  the  model  is  independent  of  wave  and 
sediment  characteristics,  and,  because  it  is  a  statistical 
method  it  requires  a  large  data  base  for  any  site  of 
interest  used.  Dean  (1973)  developed  a  heuristic  model  for 
the  prediction  of  the  presence  or  absence  of  bars  which 
compared  favorably  with  laboratory  data.  It  does  not 
describe  beach  profile  evolution  but  does  propose  a  real¬ 
istic  mechanism  of  suspended  sediment  transport  in  the  surf 
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zone  . 

Dally  and  Dean  (1980)  developed  a  realistically 
based  suspended  sediment  transport  model  for  the  surf  zone 
and  the  nearshore  region.  The  horizontal  velocities  respon¬ 
sible  for  sediment  transport  were  found  to  be  the  first- 
order  sinusoidal  component  and  the  second-order  mean  flow 
developed  from  radiation  stress  considerations.  An  exponen¬ 
tial  concentration  profile  governed  by  turbulent  and  wave- 
induced  shear  stress  was  developed  and  the  net  suspended 
sediment  flux  was  found  by  integrating  the  product  of  the 
velocity  and  concentration  profiles  over  depth.  The  results 
indicated  reasonable  quantitative  and  good  qualitative 
agreement  with  experimental  data. 

2.3  Review  of  Laboratory  Measuring  Techniques 

A  review  of  the  laboratory  measuring  experiences  by 
former  researchers  was  found  helpful  due  to  the  extensive 
laboratory  work  in  this  study.  In  the  following  sub-sections 
the  experimental  techniques  and  apparatus  involving  water 
elevation,  fluid  velocity,  and  suspended  sediment  concen¬ 
tration  measurements  are  discussed. 

2.3.1  Water  Elevation  Measurement 

In  general,  the  surface  water  elevation  measurement 
includes  the  using  of  capacitance  and  resistance  type  wave 
probe,  pressure  sensor  and  sonar  system  wave  gauge.  The 
pressure  sensor  and  sonar  system  wave  gauge  are  widely  used 
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in  the  field  study.  They  are  not  frequently  adapted  in  the 
laboratory  except  for  large  wave  tanks  where  big  waves  can 
be  generated. 

The  capacitance  and  resistance  wave  probes  are  very 
conunonly  used  in  the  wave  flumes.  The  convenience  of  the 
probe  is  that  they  are  more  sensitive  to  small  surface 
water  disturbances!  even  in  the  capillary  wave  range.  The 
wave  probe  is  usually  a  thin  steel  rod  which  does  not 
interfere  the  flow  field  significantly. 

For  surf  zone  measurement,  some  problems  arise  in 
the  accuracy  of  output  from  these  two  types  of  wave  probes. 
The  turbulent  water  may  splash  on  the  wave  probe  causing  an 
irregular  peak  on  the  output  voltage.  Also,  the  capillary 
force  of  water  on  probe  will  result  in  a  mean  water  line 
drift.  To  avoid  these  sources  of  error,  the  photographic 
method  offers  a  better  solution,  although  it  is  more  costly 
and  time  consuming. 

2.3.2  Velocity  Measurement 

Velocity  measurement  is  more  complicated  than  the 
water  elevation  measiirement ,  especially  in  the  turbiilent 
flow.  In  general,  there  are  five  types  of  measuring  tech¬ 
niques  being  used  with  successful  results.  The  instruments 
include  propeller  current  meters, photographic  method,  hot 
film  anemometers,  ultrasonic  doppler  current  meters  and 
laser  anemometer. 
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The  propeller  current  meters  are  commonly  used  to 
measure  quasi-steady  flows.  There  are  some  differences  in 
using  them  in  the  laboratory  wave  tanks  because  of  its  low 
dynamic  response  to  high  fluctuations.  Besides,  they  are 
usually  not  very  sensitive  to  low  speed  steady  flow  as 
pointed  out  by  Goda  (1964). 

With  photographic  methods,  the  water  particle  motion 
can  be  clearly  visualized  with  tracers  such  as  neutrally 
buoyant  particles  or  hydrogen  bubble.  The  water  particle 
velocity  can  be  measured  frame  by  frame  with  the  aid  of  a 
high  speed  movie  camera  and  a  slow  motion  movie  projector. 

Although  the  newly  improved  hot  film  and  hot  wire 
anemometer  can  be  operated  in  the  water  as  well  as  in  the 
air,  it  is  still  a  sophisticated  instrument.  Its  inconvenien¬ 
ces  include  frequent  failure  of  probe  head  due  to  over.- 
heating,  zero  velocity  drift  due  to  fluid  temperature  change 
and  accumulation  of  dirty  material  on  the  sensor.  It  can  not 
detect  the  flow  direction  when  the  velocity  vector  is  in 
the  plane  perpendicular  to  the  hot  film  axis.  The  draw  backs 
also  include  the  strong  non-linearity  in  the  correlation 
between  the  flow  speed  and  the  output  voltage  when  the  flow 
speed  varies  over  a  wide  range  and  also  the  film  is  sensi¬ 
tive  to  flows  that  are  parallel  to  the  film  axis. 

Ultrasonic  doppler  current  meter  and  laser  doppler 
anemometers  are  most  recently  developed  techniques. 
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Researchers  such  as  Tsuchiya  and  Yamaguchi  (1972),  Lee  et. 
al.  (197^)  had  reported  them  as  good  instrument  for  periodic 


and  random  waves  velocity  measurement.  However,  they  are 
still  in  a  developing  stage. 

2.3.3  Suspended  Sediment  Measurement 

Earlier  researchers  in  the  field  of  suspended  sedi¬ 
ment  study  were  handicapped  by  the  lack  of  suitable  measuring 
equipment.  Shinohara,  Tsubaki  and  Yoshitaki  (1958)  did  some 
pioneering  work  in  this  field.  Using  a  syphon  type  sediment 
sampler  in  a  laboratory  tank,  they  obtained  a  series  of 
measurements  of  vertical  distribution  of  the  mean  sediment 
concentration . 

A  very  useful  and  simple  technique  for  measuring 
the  suspended  sediment  concentration  in  water  was  introduced 
by  Hon-ma  and  Horikawa  ( 1963) .  By  utilizing  a  light  source 
and  light  sensor  placed  at  the  two  sides  of  a  10-cm  wide 
glass  wall  wave  tank,  they  measured  the  mean  concentration 
and  time-variztion  of  concentration  in  a  wave  period. 

Hattori ( 1969)  used  an  electrical  resistance  type 
instrument  to  count  the  number  of  sediment  particles  passing 
through  a  given  section  and  from  this  he  calculated  the 
suspended  sediment  concentrations.  Bhattacharya  (1971)  used 
an  electro-optical  instrument  for  in-situ  measurements  of 
sediment  in  suspension.  The  instrument  was  developed  at  the 
Iowa  Institute  of  Hydraulic  Research  and  is  called  Iowa 


Sediment  Concentration  Measuring  System  (  ISCMS  ) . 


CHAPTER  3 

THEORETICAL  FORMULATION 

Two  basic  surf  zone  properties,  the  surf  zone 
parameter  and  the  breaking  wave  similarity,  were  developed 
theoretically  in  this  chapter.  The  derivation  of  empirical 
formula  and  their  physical  implications  were  discussed  in 
detail.  Based  on  the  theoretical  background,  a  transport 
equation  was  developed  for  use  in  the  on/offshore  sediment 
transport  model.  Details  of  numerical  modeling  technique 
was  also  discussed. 

3*1  Surf  Zone  Parameter 

Consider  a  single  long-crested  wave  breaking  on  a 
beach  of  constant  slope,  say,  ,  the  cyclic  breaking 

process  can  be  approximately  divided  into  three  stages  as 
illustrated  in  Figure  3-1*  The  developing  stage  include  the 
period  when  the  wave  begins  incipient  breaking  at  the  crest 
until  it  collapses  on  the  beach  in  the  vicinity  of  the  water 
level  unit.  During  this  stage,  the  shape  of  the  wave  evolves 
from  highly  asymmetric  curled  form  to  a  triangular  form 
with  numerous  white  caps.  Its  potential  energy  is  converted 
into  kinetic  energy  which  in  term  is  partially  dissipated 
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Figure  3-1  Illustration  of  Breaking  Process 
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as  the  remaining  portion  preserved  to  stage  the  uprush. 


In  the  second  stage,  the  flow  surges  up  to  its  runup 
limit  similar  to  a  solid  element  moving  upslope  while  the 
kinetic  energy  is  being  converted  back  to  potential  energy. 
In  the  last  stage,  the  flow  washes  down  as  a  thin  sheet  to 
the  breaking  point  to  complete  the  cycle. 

The  time  to  complete  each  respective  stage  can  be 
estimated  based  upon  a  simple  energy  model  shown  in 
Figure  3-2.  This  energy  model  is  based  on  the  premise  that 
the  fluid  element  at  the  crest  of  the  wave  during  breaking 
remains  at  the  crest  and  is  the  same  element  that  reaches 
the  runup  limit.  This  condition  is  clearly  observed  from 
recorded  slow-met ion  films. 

If  the  kinetic  energy  of  the  fluid  element  at  the 
wave  crest  at  the  breaking  point  is  assumed  to  be  propor¬ 
tional  to  the  breaking  height,  say  M.HI,  ,  then  the 
kinetic  energy  in  the  developing  stage  can  be  expressed  as 

~  «-/) 

where  Sz  and  5/  are  the  slopes  of  the  total  energy  line  and 
potential  energy  line,  respectively.  The  velocity  of  the 
fluid  element  is 
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where  St^Si  the  corresponding  time  required 

for  the  element  traveling  from  the  breaking  point  to  the 
vfater  level  limit  is 
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where  and  yC  are  constants  of  proportionality  as  defined 
in  Figure  3-2. 

During  the  uprush  and  downwash  stages,  the  crest 
element  can  be  treated  as  a  solid  body  moving  on  a  sloped 
plan.  In  such  cases,  we  have  the  following  relationship  on 


the  basis  of  Newton's  second  law: 

?  A/^  3it>ct  -  ffA.  »yV 
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where  "f  is  a.  frictional  coefficient  and  is  the  contact 
surface  area}  the  sign  corresponds  to  the  uprush 
condition  and  the  "+"  sign  the  downwash  case.  If  we  let 
the  above  equation  simplifies  to 


^  ^  y/w  dv 

+  f  '  w  sii  dt 


where  da  is  the  thickness  of  the  runup  water  sheet. 

The  free  swash  period  in  the  region  between  the 
still  water  limit  and  the  runup  limit  can  be  obtained  by 
integrating  Equation  3-5  without  the  friction  term,  ff-^  . 


which  yields, 


I 


where  is  the  initial  uprushing  velocity  at  the  still 
water  limit.  By  letting  "id /If  ,  Equation  (3-5) 

can  be  non-dimensionalized  as 


■tv.' 

Equation  (3-7)  can  be  integrated  to  obtain  the 
uprush  and  downwash  periods  which  are,  respectively, 

£xpC^  ' 

where  we  assume  S’iZroi  for  mild  beach  slope  and  ^  is 

a  runup  coefficient  as  shown  in  Figure  3-2. 

Based  upon  Equations  (3-3) »  (3-8) i  (3-9)*  the  swash 
period  of  a  single  wave  can  be  determined: 

r.-K  fn"^ 

Where  A/ 

with 
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Therefore,  the  natural  swash  period  is  found  to  be  a  func¬ 
tion  of 

To  obtain  the  numerical  results  of  natinral  swash 


period,  the  friction  coefficient  and  the  proportional 
constants  ^  were  evaluated.  The  friction  coeffi¬ 

cient,  .  by  definition,  is  a  function  of  the  bottom 
frictional  factor,  the  initial  uprushing  velocity,  the 
fluid  element  thickness  and  the  beach  slope.  For  a  typical 
breaking  wave,  ■  8.0  fps,  tj^  ■  0.3  ft,  8ino(»0.05*  and 
assuming  f  «  O.06,  the  typical  friction  coefficient  is  2.0. 
To  investigate  the  effect  of  bottom  friction.  Equation  (3-7) 
was  integrated  in  the  swashing  region  from  Sea  Water  Level 
to  Uprush  Limit.  The  equations  obtained  for  analysis  aret 


Uprush 

H'* -f  )  €XP I- 4 K'J  (Stij 
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Downvrash 

Ctnsi 

r*  =  '^'C  I-  }  (*'*'*'> 

The  non-dimensional  period*  velocity*  and  distance  are 
defined  ast 

where  is  the  maximum  uprush  distance  and 

is  the  downwash  velocity  at  Sea  Water  Level.  The 
non-dimensional  swash  period*  velocity  and  distance  are 
plotted  against  a  wide  range  of  the  friction  coefficient 
as  shown  in  Figure  3-3. 

Figure  3-3  indicated  that  all  three  quantities  are 
decreasing  with  increasing  bottom  friction  effect.  Compared 
with  the  Friction- free  condition*  the  swash  period  of  the 
typical  breaking  wave  (^  *  2.0  )  is  reduced  to  83f5i  the 
downwashing  velocity  at  Sea  Water  Level  is  reduced  to  7195 
and  the  maximum  uprush  distance  is  reduced  to  6995.  The 
Eulerian  swash  velocities  both  without  friction  and  with 
P  *  2.0  are  shown  in  Figure  3-4.  The  fluctuating  velocity 
vectors  at  each  station  are  assumed  to  be  a  cosine  function. 
It  can  be  seen  that  the  downwash  velocity  is  lower  in 
strength  than  the  uprush  velocity  if  there  is  bottom 
friction.  The  Lagrangian  swash  velocities  for  /fl  »  0  to  10 


dcMiirtfii 

limit 


3-4  Eulerian  Swash  Velocity 
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are  shown  in  Figure  3-5* 

The  proportional  constants  were  evaluated  based  on 
a  typical  plunging  breaker.  The  constant  X/  is  1.64  for  the 
typical  wave  height-to-depth  ratio  •  •>  1.0,  and  the  typical 

crest  elevation  to  height  ratio  is  O.75.  The  breaking  wave 
height  to  depth  ratio  is  the  average  value  of  12  typical 
laboratory  plunging  breakers  conducted  on  sand  beach  as 
svmunarized  on  Table  3**1* 

The  proportional  constants  /k  and  ^  are  related 
with  the  following  equations 

This  equatlonMls  derived  based  on  the  uprushing  equations 
and  the  energy  model  shown  in  Figure  3-2. 

The  maximum  rimups  on  sand  beach  of  12  typical 
plunging  breakers  are  also  shown  Qn  Table  3-1.  The  results 
indicated  that  the  average  maximum  runup  to  breaking  wave 
height  ratio,^*  is  1.40.  Since  the  sand  beach  was  well 
compacted  and  fully  saturated  during  experiment,  the  effect 
of  bottom  percolation  was  neglected.  Based  on  the  definition 
of  friction  coefficient  and  the  knowledge  of  the  propor¬ 
tional  constants,  the  swash  coefficients  in  Equation  (^-11) 
were  evaluated  for  ranging  from  0.1  to  2.5i  which  is  the 
maximum  frictional  range  of  plunging  breaker. 


I 
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Figure  3-5  Lagrangian  Swash  Velocity  vs. 
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The  numerical  swash  coefficients  are  shown  on 
Table  3-2  and  plotted  against  yd  as  shown  in  Figure  3-6.  The 
results  indicate  that  with  increasing  bottom  friction  effect^ 
the  uprushing  period  decreases  and  the  downwasning  period 
increases •  while  the  breaking  period  and  the  total  natural 
swashing  period  remain  unchanged.  The  average  swashing 
constant,  K,  is  4.81.  The  numerical  results  imply  that 
the  total  natural  swash  period  of  a  plunging  breaker  is 
primarily  a  function  of  breaking  wave  height  and  beach 
slope,  if  the  bottom  percolation  effect  can  be  neglected. 


Having  evaluated  the  swash  coefficients,  we  now 


define  a  non-dimensional  surf  zone  parameter  as  the  ratio 
of  natural  swash  period  to  the  incoming  wave  period  (T)  : 


KT 


(3-li) 


This  non-dimensional  parameter  is  similar  to  Equation  (2-10)^ 
the  general  formula  of  previously  proposed  parameter,  with 


the  following  relation t 


where  H  is  the  wave  height  at  the  toe  of  sloping  beach  and 
is  ■‘he  breaking  wave  height. 


Figure  3-6  Swash  Coefficients 
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The  physical  implication  of  can  be  explained  as 

follows.  The  is  the  surf  zone  interference  index.  When  I 
w  w 

is  small,  each  individual  wave  will  complete  the  swash  cycle 
with  little  or  no  interference  from  the  successive  waves.  The 
flow  in  swash  zone  is  mainly  oscillating  and  the  breaker  is 
of  the  plunging  type.  When  increases,  the  degree  of 
interference  from  successive  waves  also  increases j  a  circu¬ 
latory  motion  is  gradually  developing  which  will  result  in 
a  return  flow  in  the  main  water  column;  the  breaking  pheno¬ 
menon  gradually  transforms  from  plunging  to  spilling.  This 
sequence  of  events  was  illustrated  in  Figure  2-1. 

The  number  of  breaking  waves  in  surf  zone  is  less 
than  unity  for  surging  or  collapsing  breakers.  When  the 
breaking  waves  increases  more  than  unity,  the  breaker  type 
changes  gradually  from  plunging  to  spilling.  Therefore,  by 
setting  N*»  »1.0,  the  transition  value  of  for  pure 

plunging  is  0.21.  The  N*  denotes  the  number  of  breaking 
waves  in  surf  zone.  The  transition  value  of  plunging 
breaker  is  equivalent  to  &  =  1.92,  based  on  Equation  (3-17) 
which  is  in  the  neighborhood  of  the  transition  value 
suggested  by  Battjes  (  Equation  2-8  ) . 

For  leas  than  0.21,  the  broken  ^»ve  is  able  to 
surge  up  the  slope  to  the  uprush  limit  and  return  as  down- 
wash  to  the  breaking  point  before  the  next  wave  breaks.  There 
is  no  wave  interference  and  energy  loss  is  primarily  due  to 


bottom  friction,  if  the  percolation  effect  is  neglegible 


For  this  value,  the  breaking  wave  in  surf  zone  is  runup 
dominant  and  a  large  portion  of  wave  energy  is  reflected 
resulting  in  a  partial  standing  wave  condition. 

With  the  value  of  increasing,  the  number  of 
breaking  waves  in  surf  zone  increases  proportionally.  The 
broken  wave  is  not  able  to  return  to  the  breaking  point 
before  the  next  wave  breaks.  The  increasing  wave  inter¬ 
ference  causes  greater  energy  gradient  throughout  surf  zone 
due  to  internal  friction  and  successive  breaking.  The  flow 
is  set-up  predominant  and  a  large  portion  of  wave  energy  is 
absorbed  in  surf  zone  with  little  reflection. 

In  addition  to  the  breaking  wave  properties  stated 
above,  the  I  plays  an  important  role  in  characterizing  the 
breaking  wave  similarity  which  will  be  discussed  in  detail 
in  the  following  section. 

3.2  Breaking  Wave  Similarity 

For  a  wave  train  of  single  frequency  shoaling  on  a 
beach,  a  more  or  less  energy  is  being  transferred  from  the 
primary  wave  component  to  its  higher  harmonics.  The  amplitude 
dispersing  among  harmonics  coupled  with  phase  lags  due  to 
differential  shoaling  of  each  harmonic  component  results 
in  a  highly  asyrnnetric  wave  form  upon  breaking. 


I 
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After  breaking t  as  the  wave  travels  further  inshore t 
the  wave  form  tends  to  stabilize  although  it  remains  asymme¬ 
trical.  The  conventional  higher  order  wave  theories  are  no 
longer  adequate  to  describe  these  asymmetrical  wave  forms. 

A  variation  would  be  the  introduction  of  phase  angles  among 
harmonic  components  such  that  the  water  surface  fluctuation 
is  expressed  ast 

where  a^^  is  the  n^^  harmonic  amplitude,  is  the  n^^  har¬ 
monic  phase  and  ^  is  the  fundamental  frequency. 

Equation  (3-18)  most  certainly  can  be  used  to 
prescribe  the  wave  form  of  any  particular  set  of  breaking 
waves.  Yet,  it  serves  little  useful  engineering  purpose  if 
the  solution  can  not  be  generalized.  If  we  assume  that 

(1)  Each  harmonic  amplitude  is  limited  by  the  local 
water  depth,  and, 

(2)  The  phase  velocity  of  each  harmonic  is  also 
depth  limited, 

then,  we  have. 


and 
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where  is  viave  setup.  Equation  (3-18)  can  now  be  expressed 
as: 


ait  ef»  (^nnt  ) 


(3-z/) 


Where  C^S  are  non-dimensional  coefficients  and  are 

constant  phase  angles.  Thus,  the  wave  profile  remains  similar 
and  its  magnitude  is  affected  only  by  local  parameter  (  ). 


In  shallow  water,  the  horizontal  fluid  particle 
velocity  can  be  described  by  Stokes’  higher  order  theory  as: 


The  depth  averaged  velocity  term  can  be  expressed  as: 


ir -f 

C,  ^ 


Recognizing  that  the  contribution  of  the  depth  averaged 
higher  order  terms  is  a  mean  return  flow  in  the  water  column 
as  shown  schematically  in  Figure  3-7,  Equation  (3-23) 

becomes: 


vay 
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The  fact  of  having  a  similarity  solution  of 
leads  to  the  following  similarity  solution  for  UCt) 


1 


t 
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Schematic  of  Mass  Transport  and  Return  Flow  in  Surf  Zone 
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where  and  ^  are  constant  coefficients. 

Since  C  can  be  approximated  by 

we  have 

cs-^p 

The  solution  for  U£t)  now  depends  only  upon  certain  local 
parameters  and  Uft,  .  The  similarity  constants 

5^/  depth  averaged  mean  return  flow  (Jit  SC 

obtained  experimentally  as  discussed  in  Chapters  ^  and  5> 

3.3  On/offshore  Sediment  Transport  Model 

The  present  work  utilized  a  suspended  sediment 
transport  model  to  simulate  the  beach  profile  evolution  and 
the  rate  of  beach  erosion.  The  horizontal  sediment  velocity, 
expressed  as  a  time  mean  term  and  a  periodic  term,  was  deve¬ 
loped  based  on  the  similarity  model  of  breaking  wave  in  surf 
zone  presented  in  Section  3*2.  The  predominant  sediment 
suspension  mechanism  was  assvimed  to  be  turbulent  induced 
shear  stress.  The  sediment  concentration  equation  used  was 
based  on  the  results  by  previous  researchers  and  was  assumed 


-59- 


to  be  composed  of  a  time-mean  term  and  a  periodic  term.  The 
net  suspended  sediment  flux  was  found  by  integrating  the 
product  of  the  suspended  sediment  horizontal  velocity  and 
concentration  over  depth  and  time-averaged  in  a  wave  period. 

Based  on  the  derived  transport  equation,  a  numerical 
model  was  developed  to  predict  the  foreshore  slope  evolution 
as  well  as  the  erosion  rate.  To  start  the  model,  the  initial 
beach  slope,  water  level,  sediment  characteristics  and 
breaking  wave  condition  would  have  to  be  furnished.  The 
suspended  sediment  concentration  coefficients  were  cali¬ 
brated  with  both  laboratory  and  field  data. 

3.3.1  Problem  Formulation 

Consider  a  vertical  section  of  water  in  the  surf 
zone  where  sediment  is  entrained  from  the  bottom  by  passing 
breaking  waves.  The  suspended  sediment  moves  with  the  fluid 
velocity  in  this  section  on  or  offshore  till  it  falls  down 
on  the  bottom.  The  net  suspended  sediment  flux  past  this 
section  can  be  fo\ind  by  integrating  the  product  of  the 
horizontal  fluid  velocity  and  the  sediment  concentration 
profile  over  depth,  which  can  be  expressed  as: 

where  d  is  the  water  depth  measured  from  sea  water  level 
to  bottom,  is  the  water  elevation  at  wave  crest. 
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u^(z,t)  is  the  horizontal  sediment  particle  velocity,  and 
C^(z,t)  is  the  suspended  sediment  concentration.  The  coor¬ 
dinate  system  and  expressions  of  water  depth  are  shown  sche¬ 
matically  in  Figure  (3-9) 

The  suspended  sediment  transport  equation  is  based 
on  the  knowledge  of  fluid  particle  velocity  and  suspended 
sediment  concentration.  Due  to  the  great  turbulence  in  surf 
zone,  the  conventional  linear  and  higher  order  wave  theories 
no  longer  apply  in  this  region,  neither  does  the  suspended 
sediment  profiles  developed  for  use  in  the  offshore  region. 
In  this  model,  both  equations  are  based  on  the  empirical 
formula  proposed  by  the  author  and  the  previous  researchers. 

In  the  following  paragraphs  the  horizontal  velocity, 
suspended  sediment  concentration,  transport  equation  and 
continuity  equation  used  in  the  transport  model  are 
discussed. 

HORIZONTAL  VELOCITY 

The  horizontal  velocity  of  the  sediment  particle, 
once  suspended  in  the  water  column  by  breaking  waves,  is 
assumed  to  be  the  same  as  the  fluid  velocity.  The  horizontal 
fluid  velocity  in  surf  zone  was  obtained  based  on  the 
similarity  model  for  breaking  waves  in  surf  zone  presented 
in  Section  3.2.  Based  on  Equation  (3-27),  the  velocity  can 
be  expressed  as: 


Figure  3-8  Schematic  of  Coordinate  System  and  Water  Depth  in 
Transport  Equation 
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SUSPENDED  SEDIMENT  CONCENTRATION 

The  suspended  sediment  concentration  at  a  vertical 
section  of  water  column  in  surf  zone  is  both  time  and  height 
dependent.  Using  the  technique  of  separation  of  variables, 
the  suspended  sediment  equation  can  be  expressed  asi 

m  C(t)  (i-^J 

where  C^(z)  is  the  vertical  profile  of  sediment  concentration 

and  C(t)  is  a  non-negative  time  function. 

The  time  function  is  dependent  on  the  rate  of 
dissipation  and  intensity  of  turbulence.  As  illustrated  by 
Dean  (1977)  in  Figure  3-9,  the  turbulence  intensity 
of  plunging  breaker  is  high  and  periodic,  while  that  of 
spilling  breaker  is  moderate  and  uniformly  distributed. 
Accordingly,  the  time  function  of  sediment  suspension  by 
breaking  waves  follows  the  same  pattern,  which  can  be 
expressed  by  a  Fourier  Series i 

where  and  b_^  are  Fourier  constants,  ^  is  the  funda- 
c  n 

mental  frequency,  is  the  harmonic  phase  angle.  Combining 


(A)  TURBULENCE  ENERGY  DENSITY,  e.l 
FOR  LARGE  DECAY  (PLUNGING)  ^ 


(3)  TURBULENCE  ENERGY  DENSITY,  e.(t) 
FOR  SMALL  DECAY  (SPILLING)  ^ 


Figure  3-9  Turbulence  Energy  Density  for  Large  and 
Small  Decay 


Equations  3-29  and  3-30  yields 


where  C  (z)  =  C  (z)‘C  ,  and  b  =  b/ /c”  .  For  plunging 
breakers,  the  sediment  suspension  is  periodic  with  the 
passing  breaking  waves  as  indicated  by  ISCMS  test  results 
aind  the  periodic  term  in  Equation  3-31  is  dominant.  As  the 
breaker  type  changes  gradually  to  spilling,  the  sediment 
suspension  becomes  more  uniformly  distributed  in  a  wave 
cycle  and  the  time  averaged  term  becomes  dominant.  The 
measurements  of  suspended  sediment  concentration  in  surf 
zone  supports  this  phenomenum  as  illustrated  in  Figure  3-10 
and  in  Appendix  D. 


The  time  averaged  vertical  sediment  concentration 
profile  is,  based  on  the  classical  approach,  a  balance 
between  the  sediment  grain  falling  through  a  section  due  to 
gravity  and  the  grains  transported  upward  due  to  turbulent 
mixing  and  the  concentration  gradient,  or, 


Cth  fj)  •  oO  = 


-  £ 


where  Li  is  the  particle  falling  velocity  and  £  is  the 
kinematic  eddy  viscosity.  For  open  channel  flow.  Rouse  (19^9) 
derived  the  concentration  profile  as: 


(p  -  Ci-ii  ) 


I 
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where  C  (x)  is  the  sediment  concentration  measured  at  a 
m 

reference  point  from  bottom,  &  is  a  proportion  constant 
controlling  the  shape  of  the  profile,  is  the  bottom 
shear  velocity  and  is  the  coordinate  of  the  reference 
point.  In  this  model,  the  reference  concentration  is  assumed 
to  be  proportional  to  the  ratio  of  bottom  shear  velocity 
and  sediment  falling  velocity,  or, 

CS,c»)  *  Ga  ■ 


(3-3t) 


where  the  sediment  falling  velocity  is  defined  as: 

where  ^  and  are  sediment  and  fluid  densities  respec¬ 
tively,  D^q  is  sediment  median  grain  size  ,  is  particle 
drag  coefficient  and  is  a  proportional  constant  control¬ 
ling  the  sediment  concentration. 


and 


For  sediment  suspended  by  breaking  waves,  both  iU 
has  to  be  modified  to  accomodate  the  nature  of 
turbulence  in  the  surf  zone.  In  general,  there  are  two 
quantities  which  determine  the  amount  of  sediment  entrained 


into  the  water  column: 

(1)  The  wave-induced  bottom  shear  stress, and 

(2)  The  turbulence  induced  shear  stress  due  to 
breaking  waves. 

In  surf  zone,  the  second  mechanism  is  dominant. 
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Based  on  the  derivation  of  Dally  and  Dean  (1980), 
the  turbulence  induced  shear  stress  velocity  can  be 
expressed  as: 

where  is  the  rate  of  energy  dissipation  per  volume. 

Considering  steady  conditions  when  averaging  over  a  wave 
period  and  utilizing  small  amplitude  shallow  water  wave 
relationship,  Dean  (1977)  found  that, 

where  ^ind  Jlf  is  the  breaking  wave  height  to  depth 

ratio.  Combining  Equations  (3-33)#  (3-3^)»  (3-36),  (3-37)* 
the  time  averaged  suspended  sediment  concentration  equation 
can  be  expressed  as: 

The  constants  were  calibrated  with  both  laboratory  and  field 
data  measured  in  surf  zone  by  numerous  investigators 
as  summarized  on  Table  3-3* 

Based  on  the  data  on  Table  3-3,  the  average  maximum 
sediment  concentration  in  surf  zone  is  approximately 


LABORATORIES  FIELDS 
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Table  3-3  Summary  of  the  maximum  suspended  sediment 
concentration  measured  in  the  fields 
and  laboratories 


Cj„(ppm) 

2yh 

m 

Reference  ! 

_  1 

Noda 

Bh 

0.07-0.1 

I4th  JCCE 

w 

Shimada  et  al. 

0.04-0.1 

21st  JCCE 

(A) 

Q 

Irie  et  al. 

1-10^ 

0.14-0.22 

22nd  JCCE 

UJ 

M 

Pch 

Fairchild 

- 

13rd  ICCE 

Tanaka 

0.08 

22nd  JCCE 

(B) 

1 

1 

Kana 

0.08-0.33 

l6th  ICCE 

Homa  et  al. 

10^-10^ 

0.014-0.05 

list  JCCE 

Noda 

5x10^-10^ 

0.03-0.1 

l4th  JCCE 

(A) 

C/3 

W 

Hosoi 

10^ 

0.017 

20th  JCCE 

PS 

o 

Nakato  et  al. 

10^ 

0.04 

ASCE  Vol. 

h 

108.  WWl 

<s 

a: 

Kuruhara  et  al 

^  10^ 

0.3 

3rd  JCCE 

n 

< 

Hosoi  et  al. 

10^ 

0.1 

20th  JCCE 

(B) 

Sawagaki 

10^ 

- i 

0.1 

■  1 

unpublished 

Notes  j 

(1)  is  maximum  suspended  sediment  concentration 

(2)  Z“  is  height  from  the  bottom  where  C  was  measured 

(3)  h  is  water  depth 

(4)  data  in  col’omn  (A)  were  measured  in  the  offshore 
region,  (3)  were  measured  in  surf  zone 

(5)  ICCE  means  Proceedings  of  the  International 
Conference  on  Coastal  Engineering 

(6)  JCCE  means  Proceedings  of  the  Japanese 
Conference  on  Coastal  Engineering 


-69- 

1x10  ppm  in  the  laboratory  and  1x10-^  ppm  in.  the  field.  The 
median  sand  sizes  observed  were  0.2  mm  and  the  average 
beach  slopes  were  l/lO  in  the  laboratory  and  l/25  in  the 
field.  The  above  average  data  were  used  to  calibrate  the 
proportional  constants  in  Equation  (3-38) •  The  constant  ^ 
was  found  to  be  approximately  1x10"^  gram/cm^  for  laboratory 
condition  and  2x10“^  gram/cm^  for  field  condition,  the 
vertical  profile  constant  Qt  was  also  calibrated  with 
laboratory  and  field  data  measured  in  surf  zone  by  Wang 
(1979) f  Kurihara  et.  al.  (1956),  Tanaka  (1975),  and  Sawaragi 
(1979).  The  calibrated  results  showed  that  2.5  for 
laboratory  condition  and  *  1.0  for  field  condition.  The 
suspended  sediment  concentration  and  profile  used  in  the 
numerical  model  is  shown  in  Figure  3-11*  The  measured  data 
are  also  shown  in  the  figure. 

TRANSPORT  EQUATION 

Combining  Equations  (3-27),  (3-28),  (3-31)  and 
time-averaging  over  one  wave  period,  the  suspended  sediment 
transport  equation  in  a  vertical  water  column  can  be 
expressed  as: 

a^cir)  i-  xlj^ 


Figure  3-11  Suspended  Sediment  Concentration  Profile  Used 
in  the  Numerical  Model 
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where  n  denotes  ydt  t  the  time-mean  qusintity. 

The  first  term  in  Equation  (3-39)  is  the  contribu¬ 
tion  of  mean  suspended  sediment  convected  by  the  second  emd 
higher  order  horizontal  velocity,  which  was  found  to  be 
always  in  offshore  direction.  The  second  term  in  the  equation 
is  the  contribution  of  periodic  sediment  suspension  convect¬ 
ed  on  or  offshore  by  the  fluctuating  first  order  horizontal 
velocity. 


For  plunging  breaker,  the  mean  return  flow  strength 
is  weak  and  the  sediment  suspension  is  in  phase  with  the 
first  order  velocity.  Therefore,  the  fluctuating  term  is 
dominant  and  transport  is  onshore.  As  the  breaker  type 
gradually  changes  to  spilling,  the  periodic  term  of 
suspended  sediment  becomes  relatively  unimportant 
while  the  time  mean  term  becomes  dominant.  Therefore,  the 
direction  of  sediment  transport  gradually  reverses  from 
onshore  to  offshore.  Since  the  present  knowledge  is  limited 
to  the  mean  suspended  sediment  condition,  the  present 
transport  equation  is  only  valid  for  spilling  breaker 
condition  with  I^>0,9,  in  which  the  periodic  term  can  be 
neglected  and  the  transport  direction  is  exclusively 
offshore. 

By  integrating  the  transport  equation  ever  water 
depth.  Equation  (3-39)  becomes t 
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CONTINUITY  EQUATION 

Asstiming  straight  and  parallel  bottom  contour  on 
beach,  a  two-dimensional  equation  that  governs  the  change 
in  water  depth  with  time  can  be  derived  based  on  the  conti¬ 
nuity  equation  for  solids  as: 

where  the  coordinate  system  is  shown  in  Figure  3-12,  ^  is 
the  sediment  density,  and  is  the  void  ratio  of  sediment. 
For  quartz  sand,  =2.65  gram/cm^  and  ^=0.45.  Combining 
Equations  (3-^0)  and  (3-41)  yields 


Equation  (3-42)  implies  that  the  change  in  water 
depth  with  time  in  surf  zone,  as  well  as  the  sediment 
transport  rate,  are  depth-controlled.  Starting  with  steeper 
initial  slope,  the  beach  erosion  begins  with  higher  rate. 

As  the  beach  slope  becomes  milder,  the  transport  rate  slows 
down  in  water  depth's  3/2  power.  Theoretically,  the  winter 
or  storm  profile  would  reach  equilibrium  as  the  beach 
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slope  in  surf  zone  becomes  flat  (  --^  =  0  ) 


As  the  beach  profile  approaches  equilibrium,  both 
and  can  be  approximated  by  small  constants  • 

Equation  (3-^2)  can  therefore  be  simplified  to 


I3r4i) 


where  A  is  the  lump-s\im  constant  defined  ast 

MrM  _ JL 


Integration  of  Equation  (3-^3)  yields 


Equation  (3-^5)  implies  that  the  near-equilibrium 
beach  profile  is  alvirays  concave  upward  with  h  being  x*s 
2/3  power.  This  was  proved  to  be  the  typical  equilibrium 
beach  profile  svimmarized  by  Dean  (1977)«  Therefore,  the 
present  model  is  consistent  with  the  equilibrium  shape. 

Based  on  conservation  of  mass,  the  eroded  sediment 
in  surf  zone  is  artificially  deposited  in  the  offshore 
region  seaward  from  breaking  point.  The  offshore  bar  is 
created  asymmetric  to  the  surf  profile  which  is  very 
similar  to  the  laborator^^  condition. 
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3  .J . 2  NUMERICAL  COMPUTATION 

Due  to  the  complexity  of  the  governing  equation*  a 
finite  difference  numerical  model  was  utilized  to  evaluate 
the  beach  evolution  and  erosion  rate  developed  in  Section 
3*3*1«  compute  the  beach  evolution*  the  surf  zone  was 
divided  into  several  small  cells  with  equal  width 
The  change  of  water  depth  was  computed  in  each  cell  every 
n  wave  periods*  or*  At^nr  .  The  bottom  profile  was 
smoothed  every  with  a  moving  average  subroutine 
expressed  as: 

in  order  to  simulate  the  natural  smoothing  mechanism  due  to 
bottom  shear  stress.  The  cell  width  used  were  2*5  cm  for 
laboratory  beach  simulation  and  25  cm  for  field  simulation. 


The  information  required  to  start  the  model  included 
the  breaking  wave  conditions  (  7  Tf  K  ,  ) , 

the  median  sediment  grain  size  (  )  *  and  the  initial 

beach  slope.  The  initial  foreshore  profile  for  each 
niimerical  run  was  set  in  a  concaved  upward  shape  with 
desired  average  initial  slope.  The  wave  set-up  prediction 
was  based  on  the  equations  derived  by  Longuet-Higgins: 


Jh 


A 


-76- 


The  total  volume  of  sand  eroded  at  time  t  was  the 
summation  of  the  volume  change  in  each  cell  within  the 
eroding  period.  The  eroded  material  was  located  offshore 
of  the  breaking  point  asymmetric  to  the  surf  zone  profile 
to  simulate  an  offshore  bar.  During  the  numerical  beach 
evolution*  the  surf  zone  profile  was  always  concave  upward. 
The  model  would  stop  automatically  when  the  foreshore  slope 
is  flat*  indicating  that  equilibrium  condition  has  reached. 
A  flow  chart  illustrating  the  numerical  procedures  is  shown 
in  Figure  3-13 • 


Figure  3-13  Flow  Chart  for  Numerical  Calculation 


CHAPTER  4 


LABORATORY  EXPERIMENTS 

Two  sets  of  laboratory  tests  were  conducted.  The 
first  set  of  experiments  was  designed  to  obtain  the  simi¬ 
larity  constants  and  the  mean  return  flow  strength  as 
developed  in  Section  3*2.  In  the  second  set  of  experiment 
the  instantaneous  suspended  sediment  concentration  in  surf 
zone  was  measured  for  various  wave  conditions. 

All  laboratory  experiments  were  conducted  in  the 
wave  tank  located  in  the  Civil  Engineering  laboratory.  The 
wave  tank  is  24  inches  wide  and  52  inches  deep  with  an 
overall  length  of  70  feet  as  shown  in  Figure  4-1.  The  por¬ 
tion  of  wave  tank  under  study  was  50  feet  in  length  from 
the  wave  generator  paddle  to  a  plywood  partition.  The  par¬ 
tition  was  used  to  provide  a  vertical  backing  against  which 
a  plywood  plane  beach  or  a  sand  beach  could  be  formed  and 
was  posioned  in  order  that  the  major  portion  of  the  profile 
could  be  viewed  through  the  glass  observation  panel.  Except 
for  the  20  feet  long  observation  panel  on  the  side  of  the 
tank  only,  the  tank  is  constructed  of  steel. 

The  Scotch  Yoke  type  wave  generator  is  powered  by  a 
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1750  rpm,  7.5  Horse  Power  electric  motor.  The  wave  period  is 
varied  by  varying  the  angular  velocity  of  the  rotating  drive 
arm.  The  wave  height  is  primarily  a  function  of  the  rotating 
drive  arm  offset  but  is  also  related  to  the  wave  period. 

The  steel  wave  paddle  travels  in  a  vertical  plane  on 
rollers  mounted  horizontally  on  top  of  the  wave  tank. 

Two  types  of  beaches  were  used  in  the  experiments. 
The  plane  beach  was  made  of  plywood  with  1  on  10  slope 

and  15  feet  horizontal  length.  The  plywood  beach  face  was 
painted  with  polyurethane  glass  coating  to  reduce  bottom 
effect.  The  sand  beach  material  used  in  the  experi¬ 

ments  was  graded  construction  sand  with  median  diameter  of 
0.18  mm.  The  grain  size  distribution  of  the  beach  material 
is  shown  in  Figure  4-2.  This  beach  material  is  predominantly 
composed  of  quartz  with  a  minor  portion  of  dark  minerals. 

The  experimental  procedures  and  the  arrangement  of 
tests  are  discussed  for  each  set  of  tests.  The  apparatus 
used  in  the  present  laboratory  study  is  discussed  in 
Section  4.1. 

4. 1  Discussion  of  Apparatus 

The  apparatus  used  in  the  present  laboratory  studies 
included  the  air  bubble  releasing  device,  the  photographic 
equipments  and  the  sediment  concentration  measuring  instru¬ 


ments 


Figure  4-2  Grain  Size  Distribution  of  Laboratory  Beach  Material 
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The  air  bubble  releasing  device  included  pitot  tubes 
and  compressed  air  releasing  system.  The  opening  tips  of  the 
pitot  tubes  ranged  from  0.3  to  0.5  mm.  Compressed  air  was 
released  from  a  high  pressure  steel  tank  through  a  rubber 
tube.  The  size  of  the  released  air  bubble  in  water  can  be 
adjtisted  to  a  visible  size  by  adjusting  the  releasing  air 
pressure.  The  sizes  of  air  bubble  in  water  ranged  from  1.0  mm 
to  1.5  mm  during  the  tests. 

The  photographic  equipments  included  a  Paillard 
Bowlex  16  mm  movie  camera  equipped  with  close-up  lenses.  The 
film  shooting  speed  ranged  from  12  to  62  frames  per  second. 
Kodak  Tri-X  reversal  I6  mm  BdM  movie  film  was  used.  The 
motion  pictures  were  analyzed  with  Model  224A  photo-optical 
Data  Analyzer.  The  data  analyzer  is  a  l6mm  movie  projector 
with  variable  film  projecting  speed  ranging  from  2  to  24 
frames  per  second.  Films  can  also  be  projected  on  the  screen 
frame  by  frame  similar  to  a  slide  projector. 

The  Iowa  Sediment  Concentration  Measuring  System  is 
a  complete  system  for  measuring  mean  and  statistical 
properties  of  particle  concentration  in  sediment  suspensions. 
The  ISCMS  probe  consists  of  a  gallium-arsenide  light  source, 
0.062  inch  long,  O.092  inch  in  diameter  and  a  silicon 
planar  diode  with  approximately  the  same  dimensions.  They 
were  mounted  on  supports  with  their  axes  co-lined  and  a 
clear  distance  of  0.125  inch  between  them.  The  emmision 
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spectrum  of  the  light  source  was  concentrated  at  0.9  microns 
where  the  light  sensor  is  much  wider,  although  it  peaks  at 
approximately  0.9  microns.  The  optimum  bias  voltage  of  the 
sensor  is  -20  to  -5  volts,  and  thus  it  has  a  very  large 
output  impedence  which  causes  signal  attenuation  by  the 
cable  connecting  the  transducer  to  the  instrument.  Different 
calibration  curves  should  be  prepared  for  each  sand  sample . 

The  light  attenuation  is  proportional  to  the  sediment  sxirface 
area  per  unit  volume  of  suspended  material.  These  apparatus 
are  shown  in  Photo  1. 

4.2  Water  Elevation  and  Horizontal  Velocity  Measurement 

The  purpose  of  this  set  of  experiments  were  to  carli- 
brate  the  similarity  proportional  constants  developed  in 
section  3*2.  A  total  of  10  runs  were  conducted  in  the  wave 
tank.  The  incident  wave  conditions  were  arranged  such  that 
the  breaker  types  gradually  changed  from  plunging  to  spill¬ 
ing.  while  the  the  flow  patterns  also  varied  from  oscilla¬ 
tory  to  circulatory.  The  wave  periods  were  adjusted  such 
that  they  evenly  distributed  in  between  1.0  second  and  2.0 
second.  The  breaking  wave  heights  were  also  tested  in  a 
wide  range  from  7*5  cm  to  l4.0  cm.  The  laboratory  wave 
conditions  for  the  10  runs  are  shown  in  Table  4-1. 

In  each  run  both  water  elevation  and  horizontal 
fluid  velocity  were  measured  at  various  stations  of  differ¬ 
ent  water  depths.  The  stations  were  pre-marked  with  yellow 


a 


Table  4-1  Lab.  Wave  Condition  for  Similarity  Study 
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tape  for  the  convenience  of  data  reducing.  The  scaling  grids 
were  also  pre-marked  on  the  glass  panel  with  10  ca  equal 
spacing.  Measurements  were  made  at  stations  190.  210  and  230. 
The  distance  between  each  station  was  20  cm.  The  water  level 
was  kept  constant  at  32  cm  level  in  the  wave  tank  and  the 
still  water  depths  at  the  three  stations  were  3*0.  5*0  and 
7.0  cm  respectively.  The  measuring  stations  and  water  depths 
are  shown  schematically  in  Figure  4-3. 

The  water  elevations  and  horizontal  fluid  velocities 
for  each  run  were  recorded  on  movie  film  at  1  to  3  stations 
depending  on  the  flow  conditions  in  surf  zone. The  measuring 
techniques  and  detailed  procedure  for  water  elevation  and 
horizontal  fluid  velocity  testing  are  discussed  in  the 
following  sub-sections. 

4.2.1  Water  Elevation  Measurement 

The  photographic  method  was  used  for  water  elevation 
measurement.  For  each  run  the  instemtaneous  water  elevation 
evolution  were  recorded  on  movie  film  at  various  stations. 

The  movie  films  were  then  projected  on  the  engineering  paper 
which  was  taped  on  a  flat  wall.  The  water  elevations  were 
marked  frame  by  frame  on  the  engineering  paper  and  later 
scaled  to  actual  length.  Water  elevations  were  always 
measured  up  to  the  solid  water  line  while  the  white  foam  on 
top  of  the  wave  crest  was  neglected. 


Measuring  Stations  and  Definition  of  Water  Level 
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Due  to  a  partial  reflection  effect  on  plane  beach* 
the  measured  water  elevations  were  averaged  over  several 
wave  cycles  to  avoid  biased  data.  Based  on  the  preliminary 
test  runs*  it  was  found  that  the  Instantaneous  water  ele¬ 
vation  measured  for  this  study  can  be  represented  by 
averaging  the  data  over  six  wave  cycles.  The  averaged  water 
elevations  were  non-dlmenalonallzed  by  the  local  water  depth 
hs(  )  as  presented  In  Appendix^.  A  total  of  20  wave 
elevations  were  measured  for  the  10  runs. 


4.2.2  Hbrlaontal  Fluid  Velocity  Measurement 

The  air  bubble  releasing  technique  was  applied  In 
this  study.  Pitot  tubes  were  placed  at  various  stations  In 
surf  zone.  The  tip  of  pitot  tube  was  placed  approximately 
0.5  mm  above  the  plywood  beach  floor  and  approximately  2.0  cm 
away  from  the  glass  observation  panel.  The  purpose  of  this 
arrangement  was  to  minimize  the  effect  of  bottom  and  side 
friction. 

Due  to  the  limited  width  of  wave  tank*  the  breaking 
waves  In  surf  zone  was  approximately  straight  and  parallel. 
Therefore,  it  Is  a  reasonable  assumption  that  both  water 
elevation  and  horizontal  fluid  velocity  are  uniformly 
distributed  across  the  beach*  or*  there  Is  no  lateral  varia¬ 
tion  of  flow  field  In  surf  zone.  Consequently*  the  films 
obtained  through  the  glass  panel  would  well  represent  the 
flow  at  this  section. 
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The  film  shooting  speed  was  adjusted  such  that 
there  would  he  22  to  28  frames  recorded  in  one  wave  cycle. 
Since  the  wave  periods  studied  ranged  from  1.0  to  2.0  second* 
the  film  shooting  speed  was  adjusted  from  l4  frames  per 
second  to  22  frames  per  second. 

The  principle  of  air-hubble  measuring  technique 
is  explained  in  the  following  paragraphs.  As  we  release  a 
tiny  air  bubble  through  the  tip  of  pitot  tube*  the  bubble 
in  fluid  would  reach  a  terminal  rising  velocity  shortly 
after  it  emerged  in  the  fluid.  The  constant  rising  speed 
of  the  air  bubble  through  water  column  of  uniform  stream 
would  form  a  straight  path  line.  If  the  air.  bubble  rising 
velocity  (  )  were  known*  the  flow  rate  of  the  uniform 

stream  can  be  obtained  by  measuring  the  slope  of  the  air 
bubble  in  fluid.  Let  m  be  the  slope  of  the  air  bubble  line 
measured  from  the  vertical  cross-sectional  line,  the  flow 
speed  would  be  mU^. 

If  the  air  bubbles  were  released  continuously*  the 
consecutive  rising  air  bubbles  would  form  a  straight  line 
in  a  uniform>rdtx^am~as  illustrated  in  Photo  2.  For  unsteady 
flow*  the  bubble  lines  would  be  sweeping  back  and  forth  in 
response  to  an  oscillating  flow  field. 

Assuming  that  the  fluid  flow  is  in  Stoke 's  range, 
that  is*  the  air  bubble  moves  as  the  same  speed  with  fluid 


opening  of  pitot  tube  :  0.5  mm 

Flow  rate  in  open  channel  ;  0.25  m/sec 

Grade  system  on  glass  panel  :  l*xl" 

Photo  2  Calibration  of  Air  Bubble  Rising  Velocity 
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partlcle,  the  pathlines  of  air  bubble  would  represent  the 
surrounding  fluid  flow. 

A  dynamic  analysis  of  air  bubble  in  unsteady  wave 
field  was  made  to  study  the  response  of  air  bubble  in 
unsteady  flow.  Details  of  this  analysis  is  presented  in 
Appendix  B.  It  was  found  that  •  based  on  linear  wave  theory* 
the  error  of  measured  fluid  velocity  using  al>r  bubble  tech¬ 
nique  is  less  than  20^  at  still  water  level  due  to  vertical 
acceleration.  It  was  also  found  that  the  released  air  bubble 
in  laboratory  condition  would  reach  terminal  velocity  in 
approximately  0.04  second*  which  is  2  to  of  the  testing 
wave  periods.  In  theory*  the  released  air  bubbles  would 
travel  less  than  10^  of  the  total  water  depth  before  reach¬ 
ing  the  terminal  velocity* 

In  surf  zone  where  shallow  water  condition  applies* 
it  is  reasonable  to  make  the  following  assumptions i 

(1)  The  horizontal  fluid  velocity  is  uniformly 
distributed  over  water  depth* 

(2)  The  vertical  component  of  fluid  velocity  can 
be  neglected  during  measurement. 

The  instantaneous  mean  horizontal  velocity  is  then  the 
product  of  the  air  bubble  rising  velocity  and  the  slope  of 
bubble- line  at  that  same  instant.  Figure  4-4  is  sample 
bubble  lines  reduced  from  a  test  run.  The  approximate 


\ 


straight  lines  substantiated  the  validity  of  air  bubble 
releasing  technique. 

To  obtain  an  instantaneous  horizontal  fluid  velocity 
the  movie  films  were  projected  frame  by  frame  on  the  engi¬ 
neering  paper  taped  on  wall.  The  air  bubble  line;  in  each 
frame  was  traced  on  the  engineering  paper  and  the  horizontal 
velocity  was  reduced  from  the  slope  of  the  bubble  line  at 
that  instant.  Velocities  for  six  wave  cycles  were  reduced 
and  averaged.  The  average  velocities  were  non-dimensional- 
ized  by  the  local  phase  velocity  )  and  presented 

in  Appendix  A.  A  total  of  20  velocity  measurements  were  made 
for  the  10  runs. 

4.3  Suspended  Sediment  Measurement 

The  instantaneous  suspended  sediment  concentrations 
in  surf  zone  were  measured  using  Iowa  Sediment  Concentration 
Measuring  System  (ISCMS).  This  instrument  was  applied 
successfully  in  wave  field  of  offshore  region.  Some  diffi¬ 
culties  were  confronted  as  it  was  applied  in  surf  zone. 

During  the  measurement  in  surf  zonet  dirt  and  other 
impurities  contained  in  beach  sand  were  stirred  in  the  fluid 
and  remained  in  suspension.  The  contaminated  fluid  intro¬ 
duced  higher  concentration  reading  than  the  actual  amount. 
This  kind  of  error  can  be  minimized  by  constant  washing  of 
beach  sand.  The  zero  reading  in  clear  water  kept  changing 
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from  time  to  time  due  to  the  rise  and  fall  of  water  temper¬ 
ature.  In  consequence,  zero  readings  had  to  be  adjusted  for 
each  set  of  measurement. 

There  were  air  bubble  intrusion  in  the  water  coliunn 
due  to  vrave  breaking.  The  response  of  air  bubbles  in  the 
probe  sensing  column  was  extra  light  reflected  on  the  light 
sensor,  causing  a  high  frequency  fluctuation  superimposed  on 
the  true  concentration  reading.  This  high  frequency  noise 
was  removed  by  using  a  built-in  analog  filter  with  time 
constant  O.058  second  (17.24  Hz). 

The  calibration  curve  of  ISCMS  output  voltage  vs. 
sediment  concentration  was  also  obtained.  During  the  cali¬ 
bration,  sand  sample  from  wave  tank  was  put  in  a  cylindrical 
pyrex  glass  container  with  8  inch  diameter  and  18  inch 
height.  An  electric  motor  driven  propeller  was  located 
approximately  1  inch  above  the  bottom  of  the  container. 

With  adjusted  motor  speed  and  various  amount  of  s£md  in  the 
container,  the  suspended  sediment  concentration  also  varied. 

The  ISCMS  voltage  reading  was  recorded  at  the  same 
time  while  suspended  sediment  sample  was  drained  through  a 
side  wall  of  container  into  the  50  ml  sampler.  Sediment 
concentration  was  obtained  by  weighing  the  dry  weight  of 
the  sample.  The  sediment  concentrations  were  plotted  against 
the  corresponding  ISCMS  voltage  reading  as  shown  in  Figure4>5. 
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The  calibration  curve  indicated  that  a  linear  relationship 
existed  when  the  suspended  sediment  concentration  was 
greater  than  20  gram/liter. 

During  wave  tank  measurement,  both  ISCMS  probe  and 
the  capacitance  wave  probe  were  located  in  parallel, 
approximately  5  nun  above  sand  bed.  Synchronized  sediment 
concentration  and  water  elevation  outputs  were  recorded  on 
a  two-channel  HP  chart  recorder.  The  incident  wave  heights, 
breaking  wave  heights  and  wave  periods  were  also  recorded. 

A  total  of  four  runs  were  conducted.  The  arrangement  of 
wave  conditions  are  summarized  on  Table  4-2.  The  measured 
results  are  presented  in  Appendix  C. 
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CHAPTER  5 
RESULTS 

In  this  chapter  the  data  analysis  procedures  and 
the  results  are  presented.  The  calibrated  empirical  formula 
are  applied  in  the  on/offshore  sediment  transport  model  to 
predict  the  beach  profile  evolution.  Section  5*1  presents 
the  calibrated  results  of  breaking  wave  similarity  model. 

The  experimental  data  of  previous  researchers  are  analyzed 
using  the  same  procedures  and  compared  with  the  present 
results.  The  physical  implication  and  the  region  of  validity 
of  the  similarity  model  are  discussed.  In  Section  5*2  the 
suspended  sediment  concentration  data  are  analyzed  using 
spectral  analysis  and  data  averaging  techniques.  Physical 
importance  of  the  analyzed  data  is  discussed  while  no 
quantitative  conclusions  are  made.  Section  5*3  presents  the 
niimerical  results  of  on/offshore  model  developed  in  Section 
3.3*  Both  beach  profile  evolutions  and  erosion  rates  are 
predicted  for  various  wave  and  beach  conditions.  The 
numerical  results  are  compared  with  the  laboratory  results 
of  both  the  author  and  the  previous  researchers. 

5. 1  Breaking  Wave  Similarity 
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The  20  sets  of  laboratory  water  elevation  and  hori¬ 
zontal  fluid  velocity  data  measured  as  shown  in  Appendix  A 
are  analyzed.  Detailed  data  analyzing  procedures  and  results 
are  discussed  in  the  following  paragraphs. 

In  the  first  step,  the  local  water  depths,  including 
sea  water  level  and  wave  set-up, were  reduced  from  the  movie 
film  together  with  the  water  elevation  measiarement .  The 
depths  obtained  at  each  station  were  used  to  non-dimension- 
alize  the  measured  water  elevation  and  velocity  data.  Due  to 
the  importance  of  total  water  depth  for  similarity  model, 
the  measured  set-ups  at  different  wave  conditions  were 
compared  with  the  theoretical  resiilts  by  Longuet -Higgins 
and  Stewart  (1963).  Wave  set-up  values  from  Runs  5*6  and  7 
were  compared  with  the  theoretical  predictions  as  shown  in 
Figures  D-1,  D-2  and  D-3  in  Appendix  D.  The  definition 
sketch,  experimental  results  and  the  equations  used  for  wave 
set-up  prediction  are  also  included  in  Appendix  D.  Since  the 
measured  set-ups  were  in  good  agreement  with  Longuet-Higgins 
prediction,  it  was  applied  in  the  on/offshore  transport 
model. 


For  similarity  analysis,  the  harmonic  amplitudes 
and  phase  angles  were  obtained  through  Fourier  Series  analy' 
sis.  The  principle  of  Fourier  Analysis  is  explained  briefly 
as  follows.  In  theory,  any  periodic  function  can  be 
expressed  by  the  Fourier  Series  defined  asi 


<»•+  S§  (ancos-^t  *k, 

where  f(t)  is  the  periodic  function,  t  is  the  independent 
variable  indicating  time,  T  is  the,  1st  harmonic  period  and 
®0’  ^n'  ^n  Fourier  Coefficients  defined  as  follows* 

«*» 

a-  -fa) 

These  Fourier  Coefficients  can  also  be  expressed  in 
numerical  form  as* 

a,  -  -k-'z 

n»o 

■/  Cn•At^^  cos  ±)‘Ot 
A.-  fCr,.i.-fiSi»(^jt)-ot 

nso  ' 

where  and  M  is  the  data  points  in  one  wave  cycle. 

In  this  study  the  data  points  ranged  from  22  to  2?  depending 
upon  the  wave  period  and  film  speed. 

Based  on  Equation  (3"21),  the  non-dimensional  water 
elevation  can  be  expressed  as: 

cs-z) 


i 
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The  non-dimensional  water  elevations  can  also  be  expressed 
in  Fourier  Series  aai 

(thi) 

Since  the  mean  water  level  of  each  set  of  data  was  set  at 
d+^  «or,  MWL«0,  the  time-mean  Fourier  Coefficient  is  zero 
(  Sq  *  0  ) .  To  obtain  the  Fourier  Coefficients  of  the  non- 
dimensional  water  elevation  and  velocity  functions*  data 
points  of  these  functions  were  digitized  in  a  constant  time 
interval.  The  data  points  were  then  Fourier  Transformed  with 
computer  program  for  desired  Fourier  Coefficients.  The 
amplitude  coefficients  and  phase  angles  were  determined 
based  on  the  trapezoid  rule  and  Equations  (5-2)  and  (5-3)* 


/  Af  bn  ‘ 

^  ^  An 

The  Fourier  analysis  of  horizontal  velocity  was 
based  on  Equation  (3-27) > 

Before  applying  Fourier  Analysis*  the  proportion  constant 
in  Equation  (5-6)  was  determined.  Physically*  the  value 
is  the  phase  velocity  of  breaking  wave  in  surf 
.  The  constant  'H  varies  with  different  types  of  breakers. 


zone 


for  plunging  breaker  without  interference  of  breaking  waves 
in  surf  zone  and  the  ratio  of  H^d  is  high,  a  high  "A,  value 
is  expected.  On  the  other  hand,  a  low  it  value  is  expected 
for  spilling  breaker,  since  breaking  waves  are  interacting 
with  each  other  and  the  ratio  of  H^/d  is  low. 

In  laboratory  measurement,  the  phase  velocity  of  the 
breaking  wave  is  the  siim  of  maximum  onshore  horizontal 
velocity  and  the  mean  return  flow  speed.  The  mean  return 
flow  is  just  the  numerical  average  of  the  horizontal  velo¬ 
city  over  one  wave  period,  which  is  always  in  the  offshore 
direction.  The  experimental  results  of  the  mean  return  flow 
strength  and  the  constant  ^  are  shown  in  Figures  5-1  and 
5-2  respectively.  It  appears  that  the  mean  return  flow 
strength  increases  with  increasing  I„  but  gradually  reaches 
a  constant  in  the  spilling  breaker  region.  The  ^  value,  on 
the  other  hand,  decreases  with  increasing  I  .  Since  is 
actually  the  ratio  of  wave  celerity  and  ,  the 

results  are  the  consequence  that  spilling  breakers  travel 
slower  than  plunging  breakers. 


The  wave  amplitude  coefficients,  and  wave  phase 
coefficients,  ,  of  the  first  five  harmonics  obtained  for 
the  10  runs  are  shown  on  Table  5-1*  These  values  were  plot¬ 
ted  against  the  surf  zone  parameter  as  shown  in  Figure  5-3* 
The  phase  angles  were  adjusted  so  that  ^-0.  These  data  are 
quite  scattered  but  the  trend  is  clear.  Both  and  /  become 
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Fi^re  5-1  Return  Plow  vs. 


Table  5-1  Wave  Harmonic  Coefficient 


Figure  5-3  Wave  Amplitude  Coefficients  and  Phase 
Coefficients  vs.  I 

W 


constant  when  becomes  large.  Since  a  large  corresponds 
to  a  spilling  breaker*  the  resitlts  seem  to  indicate  a  simi¬ 
larity  solutiol^  if  exists*  would  be  in  the  spilling  breaker 
region  where  is  larger  than  0.9. 

In  the  case  of  horizontal  velocity,  the  amplitude 
coefficients*  fin  .  and  phase  coefficients*  ^  *  are  shown 
on  Table  5-2  ^d  plotted  against  the  parameter,!  ,  shown  in 

W 

Pigvire  5-^»  The  phase  angles  of  horizontal  velocity  were 
also  adjusted  so  that'^=0.  These  data  reveal  that  the 
similarity  solution  also  yields  reasonable  approximations  for 
large  I^.  In  fact*  the  similarity  approximations  for  velocity 
seems  to  extend  to  regions  of  smaller  I^  than  the  corres¬ 
ponding  region  where  the  similarity  solution  of  ^  becomes 
valid*  i.e.*  velocity  profile  approaches  similar  form  as 
early  as  the  transition  flow  zone*  whereas  the  surface 
profile  will  not  become  similar  until  spilling  breaker 
condition  is  attained. 

The  experimental  results  of  other  researchers  are 
also  analyzed  with  the  same  technique  and  compared  with  the 
present  data.  Figure  5-5  plots  the  non-dimensional  surface 
variation  for  various  experients.  In  there,  data  from 
Svendsen  et.al.  (ISVA*1978)  and  Sakai  and  Iwagaki  (1978)  are 
also  shown.  Both  Svendsen *s  and  Sakai's  data  are  well  within 
the  spilling  breaker  region  (I  remging  from  1.12  to  2.27). 

W 

Their  test  conditions  are  listed  in  Table  5-3. 
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Table  5-2  Velocity  Harmonic  Coefficient 


Amplitude 


.12  .04  .04 

.10  .06  .04 


Run 

No. 

Station  . 

PB-1 

190 

210 

PB-2 

190 

210 

PB-3 

190 

210 

PB-4 

190 

210 

PB-5 

190 

210 

230 

PB-6 

190 

210 

230 

PB-7 

190 

210 

230 

PB-8 

190 

PB-9 

190 

PB-10 

190 

.01  .05 

.03  .01 


.01  .04 

.04  .02 


.06  .02 

.09  .05 

.06  .04 


.004  .02 
.02 
.03 


.05 

.03 

.05 


.06  .02 


.02  .04 


.02  .05 


-.06  -.0? 
.10  .07 


.14 

.01 


.17 

.07 


.05  -.03 

.00  .09 


.78  .90 
.99  1.03 
.90  1.14 


.01  .21 
.00  1.11 
.17  1.13 


.01  .40 

-.13  .99 

.22  .19 


.27  .91 


.19  1.24 


-.30  .81 


Figure  5-5  Non-dimensional  Surface  Profiles 
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Table  5-3  Test  Conditions  of  Svendson's  and 
Sakai  and  Isa^ki's 


Lab.  Conditiona  fISVA) _ 12Z& 


Run  No. 

Breaking  Wave 
Height,  H^(cm) 

Wave  Period 

T  (sec) 

tanac 

070703 

4.12 

1.43 

1.50 

070705 

4.43 

1.43 

1.64 

451015 

5.25 

2.22 

1.12 

451018 

4.95 

2.22 

1.09 

Lab.  Conditions  (Sakai  and  Iwagaki)  1978 


Run  No. 

Breaking  Wave 

Wave  Period 

Height,H^(cm) 

T  (sec) 

g^T'tanoc 

2-2-2 

9.4 

1.24 

1.58 

3-2-1 

11.5 

1.24 

2.27 
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All  the  data  seta  seem  to  exhibit  a  gross  similarity 
feature*  A  closer  examination  of  these  profiles  reveals  that 
the  fine  features  of  similarity  vary  for  data  sets  obtained 
by  different  investigators.  These  variations  can  best  be 
explained  with  the  aid  of  Table  where  the  statistics 
of  the  coefficients  in  the  similarity  solution  are  tabulated. 
First  of  all«  the  standard  deviations  of  the  amplitude 
coefficients  are  all  very  small  (no  more  than  from  the 
mean) •  which  means  the  similarity  solution  is  good  in  the 
region  specified.  On  the  other  hand,  the  absolute  value  of 
are  different.  The  present  laboratory  value  and  Sakai 
and  Iwagaki's  values  are  very  close.  They  both  differ  some¬ 
what  from  ISVA's  values. 

The  relative  Importance  of  each  harmonic  component 
can  be  assessed  from  the  values  listed  under  .  For 

instance,  the  fundamental  component  has  a  value  of  0.4  to 
0.5.  The  total  contribution  due  to  harmonics  higher  than 
third  is  usually  less  than  20^.  In  terms  of  wave  height,  the 
effect  of  these  higher  harmonics  is  insignificant  because  of 
the  phase  shift.  In  terms  of  energy,  their  contributions  are 
even  less.  The  relative  importance  of  the  first  three  har¬ 
monic  components  is  quite  consistent  among  the  three  sets  of 
data. 

The  results  of  phase  angle  are  somewhat  unexpected. 
All  the  data  show  that  the  second  harmonic  leads  the  first 
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Table  5-^  Statistics  of  Wave  and  Velocity  Harmonics 


AMPLITUDE  HARMONICS 
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IB 
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■ 

ISVA 

<X2 

0.11 
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X3 
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^3 
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0 
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Amplitude  ( jS  ) 

Phase 
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Mean 

O' 
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<r 
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0.078 

0.63 

- 

- 
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f«2 

0 

• 

VJ 

0 

0.054 

0.23 

t2 
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0.09 

8 

rfl3 

0.08 

0.043 

0.06 

-tj 

o.o7Tr 

0.15 

(7-  J  standard  deviation 
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harnonic  and  the  third  harmonic  leads  both  first  and  second. 
One  would  expect  the  other  way  around  as  the  higher  harmo¬ 
nics  should  have  smaller  phase  velocities.  The  values  of 
phase  coefficients  are  very  close  for  the  Delaware  and  ISVA 
data.  They  both  differ  from  Sakai  and  Iwagaki's  value.  This 
point  is  further  examined  later. 

Figvire  5-6  presents  the  non-dimensional  velocity 
profiles.  The  statistics  of  velocity  harmonics  are  also 
tabulated  in  Table  5-4.  In  here,  the  fundamental  component 
lags  the  successive  lower  component  as  expected.  However, 
these  phase  lags  are  all  very  small.  The  combined  effects 
of  dominant  fundamental  components  and  small  phase  shift 
result  in  a  more  symmetric  profile  than  that  of  wave  forms. 
The  same  kind  of  resvilts  have  been  obtained  by  Flick  (1978) 
and  Iwagaki  (1970).  A  tentative  explanation  has  been  offered 
on  the  reasoning  that  the  bottom  and  internal  friction  tend 
to  dampen  the  higher  wave  components  at  a  faster  rate.  The 
experimental  data  also  show  that  the  peak  horizontal 
velocity  always  lags  the  peak  water  surface  variation  such 
as  illustrated  by  an  example  in  Figure  5-7* 

As  the  wave  breaks  on  a  slope,  the  motion  appears 
to  be  quite  disorderly.  When  the  breaking  wave  proceeds 
upslope,  it  gradually  regains  its  regular  appearance.  The 
question  is  how  far  shoreward  from  the  breaking  point 


can 


the 


similarity  solution  be  considered 
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adequate.  For  this  purpose,  the  phase  angles  are  plotted 
against  the  non-dimensional  depth  d/d^  with  d^^  the  depth  at 
the  breaking  point.  Four  sets  of  data  from  different  inves¬ 
tigators  are  shovm  in  Figure  5-8.  In  here,  the  phase  coeffi¬ 
cients  are  defined  in  terms  of  COSINE  phase  lag,  i.e., 

\‘Z  O(r,COS  +  )  (s.^} 

n*!  * 

The  relationship  between  ^  and  ^  ,  which  is  defined  as 
the  SINE  phase  lag,  is 

t  <9-8) 

These  phase  coefficients  are  shown  to  vary  slowly  with  d/d^ 
and  approach  limits  at  intervals  when  d/d^-^0.  Since 
the  profiles  measured  by  Sakai  and  Iwagaki  are  near  the 
breaking  point,  the  results  in  this  figure  explain  why  they 
differ  from  that  of  the  present  experiment  and  that  of  ISVA; 
both  were  carried  out  in  the  inner  surf  zone  where  6./di^Q,7> 

To  utilize  the  similarity  equations  to  describe  the 
breaking  wave  properties  one  must  have  knowledge  of  and  k 
as  appeared  in  Equation  (5-6).  At  present,  we  are  unable  to 
predict  either.  The  experimental  results  of  these  quantities 
are,  however,  shown  in  Figures  5-1  and  5-2  respectively.  It 
appears  that  the  return  flow  strength  increases  with  increa¬ 
sing  I  but  gradually  reaches  a  constant  in  the  spilling 

Yt 

breaking  region.  The  k  value,  on  the  other  hand,  decreases 
with  increasing  I  .  Since  k  is  actually  the  ratio  of  wave 

IT 


9- 

are  the  consequence  that 


spilling  breakers  travel  slower  them  plunging  breakers. 

5.2  Suspended  Sediment  Concentration 

Both  spectral  analysis  and  signal  averaging  tech* 
niques  were  used  in  analyzing  the  measured  suspended  sedi¬ 
ment  data.  The  purpose  of  this  analysis  was  to  determine 
statistically  the  correlation  of  breaking  wave  and  the 
corresponding  sediment  suspension. 

A  total  of  1024  data  points  were  digitized  from  the 
strip  charts  of  both  wave  height  and  sediment  concentration 
records  of  each  run  for  spectral  analysis.  The  time  increment 
between  data  points  was  set  at  0.04  second,  therefore,  appro¬ 
ximately  22-33  wave  cycles  were  analyzed  depending  upon  the 
wave  period  of  each  run.  Details  of  spectral  analysis 
procedures  are  well  known  and  will  not  be  discussed  here. 

The  power  spectral  of  incident  wave  height,  breaking 
wave  height  and  sediment  concentration  for  the  four  runs  are 
presented  in  Appendix  C.  Based  on  the  results  of  spectral 
analysis,  some  important  physical  implications  are  discussed. 
For  all  runs,  the  single  and  concentrated  peak  of  incident 
wave  spectrum  indicated  a  periodic,  long  crested  wave.  For 
breaking  wave  spectrum,  more  peaks  appeared  at  2f  ,  3f  ,  ..... 

w  w 

where  f  is  the  fundamental  wave  frequency.  This  result 

C 

substantiate  the  harmonic  transfer  process  of  shoaling  waves 
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as  discussed  in  section  3.2.  The  sharp  lowering  of  peak 
amplitude  with  higher  frequency  component  also 
indicated  that  the  contribution  of  higher  wave  harmonic 
decreases  quickly.  There  were  also  multiple  peak  in  the 
sediment  concentration  spectrum,  indicating  a  multiple 
suspension  process  exists  in  one  wave  cycle.  The  major 
sediment  suspension  appeared  to  be  corresponding  to  the 
fimdamental  frequency  f  and  the  secondary  suspension  corres- 

C 

ponding  to  2f  .  This  phenomenum  was  more  pronounced  for 

C 

plunging  breaker  as  indicated  by  the  results  of  Run  3*  The 
double  peak  sediment  concentration  spectra  were  also  reported 
by  previous  researchers  such  as  Brenninkmeyer  ( 1973)  and 
Thornton  (1970).  The  secondary  peak  was  reasoned  to  be  the 
backward  movement  of  the  same  sediment  cloud  due  to  the 
breaking  wave  suspension. 

The  measured  instantaneous  sediment  concentrations 
were  a  random  process.  In  order  to  obtain  a  mean  sediment 
suspension  formula,  the  signal  averaging  technique  was  used. 
The  analog  ISCMS  outputs  up  to  120  wave  cycles  were  first 
converted  into  digital  data  points  with  fixed  sampling 
interval  4t.  The  digital  data  were  then  averaged  at  each 
instant  of  one  wave  cycle.  The  sampling  interval  was  set  at 
0.01  second  and  the  number  of  digital  points  ranged  from 
125  to  190  in  one  wave  cycle  depending  upon  the  wave  period. 


The  results  of  averaged  sediment  concentrations  over 
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20,  40,  60,  80  wave  cycles  for  spilling  and  plunging  con¬ 
ditions  are  shown  in  Figures  5-9  and  5-10  respectively.  The 
averaged  curves  became  stable  after  averaging  over 
60  wave  cycles.  The  curves  stabilize  faster  in  Figure  5-10, 
indicating  a  more  regular  suspension  process  for  plunging 
breaker.  The  averaged  concentration  curves  for  60  wave  cycles 
are  presented  in  Appendix  C.  These  average  concentration 
curves  can  be  expressed  as: 

^  Z  GiSiA/Cnar^ftOi)  CS-^J 

where  °0/2  is  the  mean  concentration  ,  c^  is  the  harmonic 
concentrations  and  is  the  phase  angle. 

It  was  found  that  the  mean  concentration  term 

is  predominant  for  spilling  breaker  while  the  periodic 

concentration  terra  is  predominant  for  plunging  breaker.  The 

resvilts  of  C  for  each  run  are  compared  with  the  corres- 
raax 

ponding  breaking  wave  heights  emd  spectral  energies  as 
shown  on  Table  5-5*  Due  to  the  uncertainty  of  ISCMS  applica¬ 
tion  in  surf  zone  the  measured  amd  analyzed  data  were  not 
applied  in  the  on/offshore  sediment  model. 


A  computer  program  was  developed  to  predict  the 


sediment  transport  rate.  An  outline  of  the  computer  program 
is  presented  in  Appendix  E.  A  total  of  11  numerical  runs 
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were  conducted  as  summarized  on  Table  5-6.  Three  runs  were 
made  to  compare  with  the  laboratory  results  conducted  by  the 
authors,  Sunsunura  and  Horikawa  (197^)  and  Saville  (1957)* 

Six  runs  were  made  to  demonstrate  the  effect  of  erosion  rate 
changes  due  to  variations  of  breaking  wave  height,  sediment 
grain  size,  2md  beach  slope.  One  run  was  made  to  simulate 
the  effect  of  profile  change  due  to  Sea  Water  Level  chzmge, 
and  one  rxin  was  made  to  simulate  the  erosion  process  in 
field  condition. 

Both  predicted  beach  profile  evolution  and  transport 
rate  were  compared  with  the  author's  laboratory  data  as 
shown  in  Figures  5-11  and  5-12.  The  results  indicated  that 
for  short  term  beach  erosion,  the  transport  model  compared 
well  with  experimental  data.  Figure  5-12  showed  that  the 
numerical  prediction  was  in  good  agreement  with  measured 
beach  profile  evolution  in  surf  zone,  although  it  was  not 
successful  in  predicting  the  shape  of  offshore  bar.  Since 
the  short  term  beach  erosion  is  predominantly  suspended 
sediment  transport  within  surf  zone,  the  present  numerical 
model  is  proved  to  be  well  representing  this  process. 

Numerical  runs  were  also  made  to  simulate  the  labora¬ 
tory  conditions  of  Sunamura  and  Horikawa  ( 1974) .  The  compa¬ 
risons  of  beach  profile  evolution  and  erosion  rate  are  shown 
in  Figures  5-13  5-l4  respectively.  The  results  indicated 

that  the  present  numerical  model  is  in  good  agreement  with 


Initial  Profile  =0.2  mm 


Comparison  of  Beach  Profile  Evolution  with  Author's  Results 


Figure  5-12  Comparison  of  Erosion  Rate  with  Author’s  Results 
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Comparison  of  Beach  Profile  Evolution  with  Sunamura 
and  Horikawa's  Results 
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laboratoxy  data  within  50  hour  initial  eroding  period.  The 
numerical  profile  reached  equilibrium  by  the  end  of  this 
period  while  the  laboratory  profile  kept  eroding  up  to  l60 
hours.  The  reason  for  this  deviation  was  due  to  the  fact  that 
the  present  model  predicts  only  the  offshore  suspended  sedi¬ 
ment  transport  in  surf  zone  (i.e.  from  breaking  point  to 
shoreward  limit  of  wave  set-up)  while  the  physical  model 
further  included  the  the  svi^sh  zone  erosion  which  is  mostly 
bed  load  transport  due  to  bottom  shear  stress.  Since  an 
extreme  long  erosion  period  does  not  exist  in  reality,  the 
present  numerical  model  is  adequate  for  storm  conditions. 

The  changes  of  erosion  rate  due  to  variation  of 
breaking  wave  height,  sediment  grain  size  and  initial  beach 
slope  are  shown  in  Figures  15*16  and  17  respectively.  The 
results  indicated  that,  by  holding  the  rest  of  the  initial 
conditions  unchanged,  the  beach  erosion  rate  increased  with 
higher  breaking  wave,  finer  grain  size  and  steeper  initial 
slope.  The  final  eroded  volume  was  greater  for  higher  break¬ 
ing  wave  due  to  the  wider  surf  zone  involved.  A  change  of 
wave  period  does  not  affect  the  erosion  rate  in  the  model. 
The  response  of  profile  evolution  due  to  change  of  water 
level  is  demonstrated  in  Figiure  18,  The  results  showed  that, 
by  adding  the  still  water  level,  the  breaking  point  would  be 
shifted  further  onshore  and  more  foreshore  beach  material 


would  be  eroded 


Effect  on  Profile  Evolution  due  to  Sea  Water  Level  Variation 
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Two  numerical  runs  were  made  to  simulate  the  proto¬ 
type  wave  and  beach  conditions.  The  first  run  was  designed 
to  simulate  the  large  wave  tank  condition  conducted  by 
Seville.  Comparisons  of  beach  profile  evolution  and  erosion 
rate  are  shown  in  Figures  5-19  and  5-20  respectively.  The 
results  indicated  fair  agreement  between  the  experimental 
and  predicted  data.  The  lower  predicted  erosion  rate  is  due 
to  the  use  of  prototype  suspended  sediment  constant,  which 
is  probably  too  low  for  large  wave  tank  condition.  Figure 
5-21  demonstrates  the  numerical  erosion  rates  under  labora¬ 
tory  and  prototype  conditions,  the  erosion  rates  were 
non-dimensionalized  by  the  eroded  volume  at  equilibrium 
condition.  The  results  indicated  a  lower  prototype  erosion 
rate,  which  is  probably  due  to  less  suspended  sediment 
concentration  in  the  field  and  the  scale  effect  of  physical 
model  in  the  laboratory. 
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5-19  Comparison  of  Beach  Profile  Evolution 
with  Saville's  Results 
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Figure  5-21  Comparison  of  Erosion  Rates  under  Laboratory  and 
Prototype  Conditions 


CHAPTER  6 
CONCLUSIONS 


Based  on  a  swash  model  of  breaking  waves  on  constant 
slope,  a  surf  zone  parameter  was  derived.  This  non-dimen¬ 
sional  parameter  was  defined  as  the  ratio  of  the  natural 
swash  period  of  a  single  long-crested  wave  breaking  on  a 
constant  slope  to  the  incident  wave  period,  or, 

Physically,  this  parameter  implied  the  number  of  breaking 
waves  in  surf  zone,  and  therefore,  served  as  a  wave  inter¬ 
ference  index. 


The  surf  zone  parameter  could  be  rearranged  as  the 
ratio  of  the  square  root  of  breaking  wave  steepness  to  beach 
slope.  The  new  form  was  similar  to  the  inverse  of  the 
previously  proposed  parameter  if  the  incident  wave  height 
was  replaced  by  the  breaking  wave  height.  The  surf  zone 
parameter  derived  from  a  theoretical  basis  offered  an  expla¬ 
nation  as  to  why  this  parameter  is  important  to  classify 
the  breaking  wave  properties.  In  addition,  this  parameter 
characterized  the  flow  field  and  defined  the  similarity 
property  in  the  surf  zone.  The  derived  transition  values 
from  surging  to  plunging  and  from  plunging  to  spilling 
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breaker  types  were  found  to  be  consistent  with  Galvin's 


laboratory  results. 

Assuming  both  the  water  surface  variations  and  the 
depth-averaged  horizontal  particle  velocities  in  the  surf 
zone  could  be  expressed  by  the  summation  of  Fourier 
components  I  similarity  solutions  were  proposed  in  this 
study.  These  solutions  were  based  on  a  pair  of  rather 
restrictive  conditions,  namely,  both  the  amplitude  and  the 
phase  velocity  of  each  harmonic  component  were  depth 
limited.  The  solutions,  however,  offered  the  advantage  of 
being  simple  and  completely  defined  by  a  few  local  flow 
parameters.  Laboratory  data  from  a  number  of  investigators 
including  those  obtained  at  the  present  study  were  used  to 
teat  the  validity  of  the  proposed  solutions. 

The  results  seem  to  suggest  that  the  similarity 
solution  is  not  universally  applicable.  However,  if  the 
wave  is  of  the  spilling  type  and  is  far  inshore  in  the 
inner  breaker  zone,  the  similarity  solution  become  suitable 
to  describe  the  mean  flow  characteristics.  These  conditions 
can  be  defined  in  terms  of  the  surf  zone  parauneter  and 
the  relative  depth  parameter  d/d^.  The  region  of  suita¬ 
bility  is  found  to  be  when  I^^0.9  and  d/d^<0.6.  In  this 
region,  both  surface  profile  and  horizontal  velocity  can  be 
adequately  described  by  the  first  three  harmonics.  The 
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velocity  profile  is  foxond  to  be  less  asymmetric  than  the 
surface  profile.  For  surface  profile,  the  higher  harmonics 
lead  the  lower  ones  whereas  the  contrary  is  true  for  the 
velocity  profile. 

The  results  of  analyzed  suspended  sediment  concen¬ 
tration  data  indicate  two  physical  significances.  The  mul¬ 
tiple  peaks  in  the  sediment  concentration  spectrum  implies 
that  miiltiple  suspension  process  exists  in  one  wave  cycle. 
The  peaks  are  more  pronoxinced  dxiring  spilling  wave  condition 
resxilting  in  the  uniform  distribution  of  sediment  concen¬ 
tration  in  one  wave  cycle.  The  averaged  instantaneous 
sediment  concentration  profile  is  periodic  and  is  in  phase 
with  the  velocity  profile  for  plxmging  breaker.  The  concen¬ 
tration  profile  for  spilling  breaker  is  more  xjniformly 
distributed  with  its  peak  lags  the  velocity  peak.  The 
difference  in  sediment  suspension  pattern  is  foxand  to  be 
the  controlling  factor  of  on/offshore  transport  direction. 

As  part  of  an  attempt  to  simulate  the  on/offshore 
transport  phenomena  in  the  nearshore  region,  the  present 
nxamerical  model  is  successful  in  predicting  the  offshore 
sediment  transport  both  qualitatively  and  quantitatively. 
Given  any  initial  conditions  including  breaking  wave 
height,  sediment  grain  size  emd  beach  slope,  the  model  is 
able  to  predict  erosion  rate  and  profile  evolution  in  surf 
zone  provided  that  the  breaker  is  of  spilling  type. 
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A  suspended  sediment  transport  equation  is  applied 
in  the  model.  The  transport  equation  is  obtained  by 
coupling  the  horizontal  particle  velocity  with  the  sus¬ 
pended  sediment  concentration.  Due  to  the  limited  Icnowledge 
of  sediment  suspension  in  surf  zone,  the  present  model  is 
only  able  to  predict  the  mean  condition,  which  is  the 
offshore  transport  case.  The  model  is  highly  depth- 
controlled.  The  erosion  rate  reduces  sharply  as  the  beach 
slope  becomes  milder.  The  equilibrium  condition  would  be 
reached  as  the  beach  slope  in  surf  zone  is  flat. 

The  numerical  model  can  be  applied  to  both  labora¬ 
tory  and  field  conditions  by  using  the  proper  proportional 
constants  which  are  calibrated  with  either  laboratory  or 
field  data.  The  present  model  predictions  are  found  t.o  be 
in  good  agreement  with  all  laboratory  results  within  50 
hour  erosion  period.  Beyond  this  period,  berm  crest 
erosion  takes  place  due  to  eroded  unstable  foreshore  slope. 
This  erosion  process  should  be  included  in  future  model 
study. 


Although  the  present  laboratory  measuring  tech¬ 
niques  offer  sufficient  engineering  estimates  of  the  flow 
field  in  surf  zone,  a  more  sophisticated  velocity  measuring 
technique  is  desirable  for  future  study.  The  water  elevation 
and  velocity  measurements  should  be  conducted  for  various 
beach  slopes  and  with  both  plane  and  sand  beaches.  More 
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detailed  laboratory  investigation  of  sediment  suspension 
in  surf  zone  is  necessary  for  a  complete  on/offshore  sedi¬ 
ment  transport  model. 
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APPENDIX  A 


WATER  ELEVATION  AND  HORIZONTAL  VELOCITY  DATA 


NON-DIMENSIONAL  QUANTITY  NON-DIMENSIONAL  QUANTITY 
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Figure  A-1 


APPENDIX  B 

DYNAMICS  OF  AIR  BUBBLE  IN  UNSTEADY  WAVE  FIELD 

The  coordinate  system  for  the  equation  of  motion  of 
air  bubble  in  fluid  environment  is  defined  as  follows: 
x-y  coordinate  for  fluid  field 
coordinate  for  air  bubbles 

The  two  coordinates  are  shown  schematically  in  Figure  B-i, 
The  vertical  equation  of  motion  of  air  bubble  can 
be  expressed  as: 

+  ”/«  )  CS-0 

where  is  air  density 

fit  is  fluid  density 

is  air  bubble  volume 
Ap  is  the  project  area  of  air  bubble 
a,  is  drag  coefficient 
ntA  is  the  added  mass  of  air  bubble 
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For  a  sphere  shaped  air  bubble  with  diameter  D,  its  volume, 
project  area  amd  added  mass  can  be  expressed  as* 


Equation  (b-1)  is  a  non-linear  second  order  differential 
equation.  Assuming  that  the  flow  field  is  in  Stokes’  range, 
the  non-linear  term  in  Equation  (B-l)  can  be  simplified  to: 


a/»/»  4> 


where 


Substituting  Equation  (B-2)  in  Equation  (B-l)  and  dividing 
each  term  in  the  equation  by  V^,  Equation  (B-l)  becomes: 


fp 


Equation  (B-3)  can  be  rearranged  toi 


)  CB-4-) 
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Since  ^  can  set  =0.  Dividing  both  sides  of 

Equation  (B-4)  with  fia  yields: 

^B-s) 

For  air  bubble  releasing  from  bottom,  the  initial  condition 
can  be  expressed  ass 


For  a  sinusoidal  fluid  particle  velocity  in  the  vertical 
direction,  we  define 


jf-  cas  (S-^J 


where  Z4t  is  the  maximum  amplitude  of  vertical  velocity 
60  is  the  wave  frequency 

Integration  of  Equation  (B-5)  applying  the  initial  condition 
and  Equation  {3-6)  yields: 


’I 


=  Vk[i-Oi 

-L 


l)‘m 


C6suii 


ZAM 


where  is  the  terminal  rising  velocity  of  air  bubble 
defined  as  ^  and  a=3  -ik- 

Equation  (3-7)  can  be  simplified  to 
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where 


The  horizontal  equation  of  motion  of  air  bubble  can 
be  expressed  as: 

Equation  (B-9)  can  be  simplified  by  assuming  0  i 


Cok 


For  a  sinusoidal  horizontal  fluid  partical  velocity,  we 
define : 


j(»6/ni  Sin  CyO± 


where  Urn  is  the  maximum  amplitude  of  horizontal  velocity. 
The  initial  condition  can  be  defined  as: 


\(o)-o  at 


Integration  of  Equation  (B-iO)  applying  the  initial 
condition  and  Equation  (B-H)  yields: 

The  analytical  solutions  of  vertical  and  horizontal 
air  bubble  velocities  are  discussed  in  the  following 
paragraphs :  _ 

letting  V=  r-  ‘***^*’^ 
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The  vertical  and  horizontal  air  bubble  velocities  in 
Equations  (B.-8)  and  (B-12)  becomes: 


5  «  n.  fo<  «■«  €~‘*  *J  (0-t4.) 


For  small  air  bubble, 

»° 


o(.= 

/J  = 

r  = 
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^  JeL 

^  A 


■as  -2. 


/■+  i/'^y 


W  \Z 


-/  2AUJ 


✓  -  *v 

C  — 


4^ 


CJ 

/ 

2.t^ 
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Cen^^ 


Therefore,  for  a  small  air  bubble  rising  through  calm  water, 
the  vertical  and  horizontal  velocities  can  be  expressed  as: 

^-Wb  c*- 

ceii) 

^  =  0 

The  traveling  time  of  air  bubble  in  calm  water  as  a  function 
of  the  size  factor  ”a"  is  illustrated  graphically  in 
Figure  B-2. 
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For  the  present  laboratory  condition,  the  air  bubble  reached 
tejrminal  velocity  in  approximately  0.06  seconds. 


For  the  case  of  air  bubble  in  unsteady  wave  field, 
the  vertical  and  horizontal  velocities  of  air  bubbles  after 
they  reached  terminal  velocities  are: 


CB-17) 


The  nximerical  values  of  oi  and  ^  are  plotted  against  wave 
period  "T"  for  various  "a"  values  as  shown  in  Figure  B-3.  It 
can  be  seen  that  for  longer  wave  period  and/or  small  bubble 
size,  /  ar.d  6  S^O  •  The  results  imply  that  for 

smaller  bubble  size  and/or  longer  wave  period,  the  better 
the  fluid  particle  velocity  can  be  represented  by  air  bubble 
velocity. 


In  order  to  find  the  maximum  error  imposed  during 
the  velocity  measurement,  a  correction  factor  "CF"  was 
introduced.  Referring  to  Figure  B-4,  the  measured  mean 
horizontal  velocity  can  be  expressed  as: 

ZT  —  ^  [ZTb -hsiVln  (6-18) 
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Letting  the  maximimi  bubble  line  slope  measured  be  m,  the 
bubble  line  slope  at  any  instant  in  one  wave  cycle  can  be 
expressed  as  ^  ^ 

JL.  .  ^  s,-)f  / 

Vh  / 

Physically,  for  small  phase  error  ",  m  5^^  is  the  true 
horizontal  velocity  at  that  instant  and  is  the 

error  introduced  by  vertical  fluid  velocity.  The  maximtun 
error  imposed  is  then  as  shown  in  non-dimensional 

form.  The  correction  factor  is  defined  as 

The  correction  factor  was  plotted  against  different  water 
depths,  wave  heights  and  wave  periods  as  shown  in  Figures 
B-5  and  B-6.  The  results  indicated  that  at  wave  trough,  CF 
ranged  from  0.3  to  0.4  for  deep  water  condition,  0.2  to  0.3 
for  shallow  water  condition.  For  fluid  particle  velocity  in 
surf  zone  the  error  is  further  reduced  due  to  weak  vertical 
fluid  velocity.  Therefore,  the  maximum  error  introduced 
using  air  bubble  velocity  measuring  technique  would  be  less 
than  20?5. 
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Pigure  B-1  Coordinate  System  for  Air  Bubble 
and  Fluid  Particle 
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B-4  Correction  Factor  vs.  Total  Water  Depth 


APPENDIX  C 


WAVE  AND  SUSPENDED  SEDIMENT  SPECTRA  AND  CONCENTRATION  DATA 
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Figure  C-1  Incident  Wave  Spectrum  for 
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Figure  G-2  Breaking  Wave  Spectrum  for 


F/e  s/s 


AD-AllS  AA2  DELANARE  UNZV  NEWARK  DEFT  OF  CtVlL  ENStNEERlNS 

SURF  ZONE  properties  ANO  ON/OFFShORE  SEDIMENT  TRANSPORT. (U) 

JUN  82  W  VANS*  H  HANS  NS00i4-81-K-0297 

itMri  *<«tFTFn  ce-a2»a5 _ _  NL 


Figure  C--5  Incident  Wave  Spectrum  for 
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Figure  C-6  Breaking  Wave  Spectrum  for 
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Figure  C-7  Suspended  Sediment  Concentration  Spectrum  for  Run 
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Igure  C-9  Incident  Wave  Spectrum  for  Run 


Figure  C-10  Breaking  Wave  Spectrum  for  Run 


Figure  C-11  Suspended  Sediment  Concentration  Spectrum  for 
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Figure  C-12  Averaged  Suspended  Sediment  Concentration  for  Run 
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Figure  C-13  Incident  Wave  Spectrum  for  Run 
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Figure  C-l4  Breaking  Wave  Spectrum  for 
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Figure  C-15  Suspended  Sediment  Concentration  Spectrum  for 


APPENDIX  D 


WAVE  SET-UP  ON  PLANE  BEACH- 
THEORY  AND  LABORATORY  DATA 

A.  Set-UT>  Equations 

Based  on  Longuet-Higgins  ( 1963) ,  the  set-up  equation 
can  be  expressed  as: 

<?;!  ,  . 


where  H^  =/fh  for  spilling  breaker  and  h^d+J^  (see  Pig.D-1). 


^BRE4Wa/4-  PoiA/T 


Figure  D-1  Definition  Sketch  of  Wave  Set-up 
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At  breaking  point i  x«0,  we  have 


At  still  water  level*  z«x^»  we  have 


/•N4> 


i.  *<4  (  77^  -  <4  f-fJV/-  -fTlf A#*; 

OC,  ■/««  -He  -Pfwt  <*••> 

=  Wi  <'/- <►<) 

B.  Laboratory  Data 

Three  sets  of  laboratory  wave  set-up  data  were 
reduced  from  the  water  elevation  data*  corresponding  to 
Rvin  No.  5*6*7  in  Table  4-1.  The  wave  and  beach  conditions  of 
these  three  runa  represent  typical  spilling  breaker.  The 
measured  set-ups  are  sianmarized  in  Table  D-1  and  compared 
with  the  theoretical  predictions  shown  in  Figures  D-2  *  D-3 
and  D-4 . 


1 


tiUre^  levBi.  can) 


-193- 


HcminoMim,  can) 

Figure  D-2  Wave  Rian-up  Data  for  Run  5 


imvmL  (cm> 


-194- 


Figure  D-3  Wave  Run-up  Data  for  Run  6 
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Figure  D-4  Wave  Run-up  Data  for  Run  7 


APPENDIX  E 


NUMERICAL  MODEL  PROGRAM  LISTING 
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100 

♦RESET  FREE 

110 

♦SET 

AUTOBIND 

120 

FILE 

6<MA,<RECSIZE»22) 

130 

DIMENSION  CC20),T<20)»XP(20> 

140 

COMMON  X<2001)»Y(2001»20)»N 

150 

COMMON/B1/G.PI,PI2 

160 

C0MM0N/B2/XB  »  TB  »  XWL » A YO » XPL  f YPL 
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C0MM0N/B4/GS » PO  »  DM 

190 

C0MM0N/B5/WML  <  2001 f  20 ) »  Q ( 20 ) 

200 

C3k«« 

210 

INPUT  NUMERICAL  CONSTANTS 

220 

C*** 

230 

6=980.62 

240 

PI=3.1416 

250 

PI2=2.#PI 

260 

6S=2.65 

270 

P0=0.45 

280 

C* 

G=GRAUITATIONAL  ACCELERATION  (CM/SFC**2) 

290 

C« 

GS=SEDIMENT  GRAIN  DENSITY  (GRAM/CMJktS) 

300 

C« 

PO=SEDIMENT  POROSITY 

310 

320 

C««3|C 

INPUT  X»Y  COORDINATE  AND  INITIAL  SI.OPE 

330 

C*** 

• 

340 

DX=2.5 

350 

YMAX=60. 

360 

WRITE(6>602) 

370 

602 

F0RMAT<1X»*TYPE  1  FOR  1/10  SLP.2  FOR  1/15 

380 

1  3  FOR  1/20  SLP'> 

390 

READ<5»/)  ISLP 

400 

IF(ISLP.EO.l)  GO  TO  1 

410 

IF<ISLP.Ea.2)  GO  TO  2 

420 

IF<ISLP.EQ.3)  GO  TO  3 

430 

1 

XMAX=600. 

440 

SUL =30. 

450 

SLP=1 ./lO. 

460 

N=XMAX/DX 

470 

DO  11  I=1»N 

480 

11 

Y<I»1)=YMAX-0.8*(DX*(I-1 > ) **0.6667 

490 

GO  TO  5 

500 

2 

XMAX=700. 

510 

SUL=35. 

520 

SLP=1./15. 

530 

N=XMAX/DX 

540 

DO  12  I=1»N 

550 

12 

Y(I»1)»YMAX-0.7*(DX*<I-1>  >**.6667 

560 

GO  TO  5 

570 

3 

XMAX-800. 

580 

SUL=40. 

590 

SLP=l./20. 

600 

N=XMAX/DX 

198 


610 

DO  13  I»1»N 

620  13 

Y<I»l)»YMAX-0.6i|t<nXi|t<I-l )  ) **0.6667 

630  5 

CONTINUE 

640 

DO  14  1-1>N 

650 

X<I+t)»X<I)+DX 

660  14 

Y<I»2)»Y<I»1) 

670  C* 

YMAX«MAXIMUH  BEACH  ELEUATION  <CM) 

680  C« 

SWL-SEA  WATER  LEVEL  <CH) 

690  C;x 

N-NUHBER  OF  DISTANCE  INCREMENTS 

700  C*** 
710  C««« 

INPUT  WATER  LEVEL »SEDIMENT  SI7E»  EROSION 

TIME 

720  C*** 
730  15 

WRITE(6»604) 

740  604 

FORMATdX,  'RN*  •»  'DM*  *.  ‘DTf  •»  'NDT.  ■  > 

750 

READ(Sf/)RN»DM»DT»NDT 

760  C« 

RN-RUN  NUMBER.  IF  SUL  CHANGED.  RN  ALSO  CHANGE 

770  C* 

DM-MEDIAN  GRAIN  SIZE  <MM) 

780  C* 

DT»DURATION  OF  TIME  FOR  BEACH  EVO.  (HOUR) 

790  C* 

NDT»NUMBER  OF  DT  DESIRED 

800  C««« 
810  C«»» 

INPUT  INCIPIENT  WAVE  INFORMATION 

820  C««« 
830 

URITE<6>606> 

840  606 

FORMAT< 1X» 'HB. • » "DB. • » 'PER. • ) 

850 

REAn<5»/)  HB.DB.PER 

860 

CKAP»HB/DB 

870  C* 

HB*BREAKING  WAVE  HEIGHT  AT  BREAKING  POINT 

(CM) 

880  C» 

DB»WATER  DEPTH  AT  BREAKING  POINT  (CM) 

890  C» 

PER»INCIDENT  WAVE  PERIOD  (SECOND) 

900 

910  C**» 

COMPUTE  BEACH  EVOLUTION  IN  nT*NnT 

920  C*»* 
930 

YBl»SWL-nB 

940 

CALL  XYCORDd.XB.YBl.IB.l) 

950 

DO  20  J*3fNDT+2 

960  C» 

COMPUTE  WAVE  SET-UP 

970 

J=J-1 

980 

CALL  WAVSUP(XBrAYB.AYO.XPL..I) 

990 

XPL*XB  -XPL 

1000 

XP< J)=XPL 

1010 

YB=YBlfAYB 

1020  C* 

FIND  N'TH  INCREMENT  AT  SEA  WATER  LEVEL 

1030 

CALL  XYCORDd.XWL.SWL.IRf J) 

1040 

IWL»=IR 

1050  C* 

FIND  N'TH  INCREMENT  AT  RUN-UP  LIMIT  AND 

YPL 

1060 

CALL  XYCORD(O.XPL.YIRf IR» J) 

1070 

YPL®YIR 

1080 

IPL=IR 

1090 

SUP»  < YPL  - YB ) / ( DX* ( I B~ I PL ) > 

1100 

JaJ+1 

1110 

DO  25  lalPL.IB 

199- 


1120  25  WML<I»J)»YPL-SUP*nX*(I-IPL) 

1130  DO  30  I=IB+1»N 

1140  30  UML<I>J)»SUL 

1150  C*  COMPUTE  EROSION  RATE  IN  DT 

1160  CALL  BEHEV0<J»N»DT»DX»PER1 

1170  T<1)=0. 

1180  T< J)«< J-2)*nT 

1190  20  CONTINUE 

1200  C*  OUTPUT  TO  PRINTER  OR  PLOTTER 

1210  URITE<6>608> 

1220  608  FORMAT < IX » 'TYPE  0  FOR  PRINT  OUT,  I  FOR  PLOTTER") 

1230  READ<5,/)ID 

1240  IFdD.EO.l)  GO  TO  40 

1250  WRITE<6»6010)  RN 

1260  6010  F0RMAT<10X»*RUN  NUMBER=",I3> 

1270  WRITE<6»6020)  HB 

1280  6020  FORMATdOX,  "BREAKING  WAUE  HEIGHT  f  CM)*"  ,F5. 1 ) 

1290  URITE<6»6030)  PER 

1300  6030  F0RMAT<10X»"MAVE  PERIOD  < SEC ) = " « F5 . 1 ) 

1310  WRITE<6,6040)  SLP 

1320  6040  F0RMAT<10X, "INITIAL  SLOPE*" ,F5,2) 

1330  WRITE <6, 6050)  DM 

1340  6050  FORMATdOX,  "MEDIAN  GRAIN  ST7F  (MM^*".F5.2) 

1350  URITE(6,6060) 

1360  6060  FORMATdOX,  "ERODING  TIMF(HOUR)  "  ,2X,  "ERODED  UOI.(Cm/j 
1370  WRITF(6,6070)  < (T< J) ,0( J> ) ,J=3,NDT+2) 

1380  6070  F0RMAT<2F20.2) 

1390  URITE<6,/)  <XP<J )  ».J=*2,NnT+7) 

1400  DO  50  J*2,NDT+2 

1410  WRITE<6,6080)  T<J) 

1420  6080  FORMAT<//, IX, "ERODING  TIME* " , 15 , // ) 

1430  URITE<6,6035)  XP<J-1) 

1440  6085  FORMATdX,  "ERODING  DTST.*  ".riO.J, //> 

1450  DO  50  K*IPL,IB,10 

1460  URITE(6,6090)  < Xd  )  *  I»K » K+9 ) 

1470  6090  FORMATdX, "X*  ",10F7.1) 

1480  WRITE<6,6100)  < Y ( I , J > , I«K ,K+9 ) 

1490  6100  FORMATdX,  "Y*  ",1,0F7.1,/) 

1500  50  CONTINUE 

1510  40  CONTINUE 

1520  C*  OPTION  FOR  SUL  CHANGE 

1530  URITE<6,610) 

1540  610  FORMATdX,  "TYPE  0  FOR  STOP,  t  FOR  SUI.  CHANGF") 

1550  RFADf5,/)  ISU 

1560  IE < ISU. EH. 0)  GO  TO  60 

1570  DO  70  1*1, N 

1580  Yd,t)»Y<T,.l> 

1590  Y<I,2)-Yd,J) 

1600  70  CONTINUE 

1610  GO  TO  15 

1620  60  STOP 


1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1360 

1370 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 


END 

SIJBROUTTNfT  WAU8MP^VB,AYR.AY0,XPI.  .  J) 

ClicKc* 

C*  THIS  SUOPOlITTNF  COMPOTFS  UAUF  srTIJP  IN  SIJRF 
C*  AYB=»SFTnp  at  HRFAKINH  POINT  ^ABOUF  BUI  > 

C*  AYO=SET|JP  AT  SUI  (ABOUF  BUI  1 
C*  XPL.=X-CnRri.  AT  RUNUP  I  TMTT 
COMMON  X<2001 1 .y<2001 .201 
COMMON/BB/TPI  .TB.TUI  .BUI  .CKAP 
CAI.  I  XYCORn<O.XB.YTR.  TR..n 
riB=BUI.-YTR 

CAI  I  XYCORBM  .XTRN.BUI  ,TR..n 
S»nB/(XB-XIRNl 

AYB*B*(XB-XTRN1/I1  . +CK:AP**2/1  6  -  ) 

AYOsAYR*^  t  .-5.)ltCK  AP*!fc2/1/..  1*^1  .  I  rkAP)k*2/1  6  *  > 

1  *  ^  B  .  UcCKAPiktl?/ 1  6  .  > 

XPI.=  (XB-XIRNlilt^  1  .+FXP/-0.6*M- 21  1  *?%  .  *CI<AP**2/ 1 , 

RETURN 

ENB 

C*** 

C*#* 

SUBROUTINE  BFHFUO  ^ .  I . N . BT . BX . PFR ) 

C!k>kilc 

C*** 

C*  THIS  SUBROUTINE  COMPUTES  TRANSPORT  RATES 

C*  AND  BEACH  FUOI..UTTON 

DIMENSION  D0(2001.201 .nY(2001 .201 
COMMON  X( 2001 >«Yr 200 1.201 
C0MM0N/B.1/0.PT.PT2 
C0MM0N/B2/XB.YR.XUI  .AYO.XPI  .YPI 
C0MM0N/.B3/TPI  .TR.TUI  .SUI  .CKAP 
COMMON /B4 /OS . PO . DM 
C0MM0N/B5/UMI  ( 2001 . 20 1 . 0 C 20 1 
C*  INPUT  NUMBER  OF  UAUE  PERIODS  EACH  ITERATION 

URITE(6./>12) 

612  FORMAT M X. ■NO.  OF  UAUF  PERTOBS  EACH  TTER.’l 
RFAB(5,/1  CIN 
NIT=TNT(3A0.*BT/<PEP*CTN1 1 
BO  10  1=1.1  .N 
10  Y(T.  J1=Y(T..I~1  1 

BO  15  T  =  1.TPI-1 
15  UMI  ^T..I)=YfT..ll 

C*  COMPUTE  TRANSPORT  EOllATTON 

AF*1  . 
r.K=r0.01 
Cll«0.25 
CBal .5 

UF=SnRT<4,*R*<BM/10. 1*(RS-1 . 1 /(S.tCBl 1 
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2140 

2ir^o 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 


ACK  ( 1=  (  PK  *r  KAP*r  1 1*04*  1  .  5  )  /  M  .  *AF*UF**2  ) 

nn  20  K=1»NIT 

no  20  T=TPI  +1 .TR-1 

onwi  Y(  T-1  ,.n-Y(  T+1  ..n  )  /  < 2**ny  ^ 

TF< AOHl .1  F.o. >  AnHIaO.OOOOl 

AnH2=a^  fUMI  ^T  +  1  ..0~Y^T+1  ..1)  )44t  .5-fU«l  <'T-1  ..I) 

1  -Y(  T-1  ..m**!  .5> 

AFXP=FXPr-2.!ltAF*OF/RORT(CKAP*r5!(CAnH1  1  1 
nn^  T  •.Jl-0.54Ar,KII«M  .-AFXPIA^  AnH2*AnH1  rTN*PFR) 
Q<  J)=0<.ll+no<T..Jl 
OY  <  T » .  1 1 =00  ^  T » .  n  / ^  nx*RR4 ( 1  * -pn  ^  s 
Y  f  T . .  n  =Y / T . J ) -OY  f  T . .  n 
Yr  ip».  »,n=Yr  TPi  +1  ,.n 

n*  SMOOTH  THF  OFACH  PPOFTI  F  IISTNR  HANNTNO  UTNOOU 
CAI  .I..  SMOOTH  Mi 
20  CONTTNIJF 

C*  SET  UP  OFFSHORE  RAR  ASYMMFTRTr  TO  RFAPH  PROFT!  F 
IPPaTR-TPI  -1 
DO  30  T=1.TRP 
OnY=Y^ TR-T .1-1 >-Y^  TP-T.  n 
Y(  TB+I»  J)=Y^  TR+T..I)+nOY 
30  CONTINUE 

RETURN 
END 

C444 

C#** 

SUBROUTINE  X YCORD < IND . XTR . YTR . TR , , n 

C444 

C444 

C4  THIS  SUBROUTINE  FINOS  Y-CORO  RTUFN  X-CORO 
C*  AND  VICE  UFRSA 

C*  IF  INn=0»SiyFN  X-CORO  TO  FINO  Y-CORO 

C*  IF  TN0=1.GIUFN  Y-CORO  TO  FIND  X-CORO 

COMMON  X<2001  i  .YMOOl  .20i  ,N 
IF<TN0.E0.1)  no  TO  30 
DO  10  1=1 .N 

IF< <XTR-X( li i .LE.O. )  00  TO  20 
10  CONTINUE 

20  TR=T 

YIR=Yf  IR. ..n 
RETURN 

30  DO  40  1=1 »N 

IE< ( YTR-Y< I . J) ) .OE.O. >  00  TO  50 
40  CONTINUE 

50  IR=I 

XIR=X<TRi 

RETURN 

END 

C«»« 

C444 


SUBROUTINE  SMOOTH <.!> 
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